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^ ■ ABSTRACT 

p • We study Big Crunch/Big Bang cosmologies that correspond to exact world- 

ly , sheet superconformal field theories of type II strings. The string theory space- 

time contains a Big Crunch and a Big Bang cosmology, as well as additional 
"whisker" asymptotic and intermediate regions. Within the context of free 
^ I string theory, we compute, unambiguously, the scalar fiuctuation spectrum 

■ in all regions of spacetime. Generically, the Big Crunch fiuctuation spectrum 
is altered while passing through the bounce singularity. The change in the 
spectrum is characterized by a function A, which is momentum and time- 
dependent. We compute A explicitly and demonstrate that it arises from 
the whisker regions. The whiskers are also shown to lead to "entanglement" 
entropy in the Big Bang region. Finally, in the Milne orbifold limit of our 
superconformal vacua, we show that A — 1 and, hence, the fiuctuation spec- 
trum is unaltered by the Big Crunch/Big Bang singularity. We comment on, 
but do not attempt to resolve, subtleties related to gravitational backreaction 
and light winding modes when interactions are taken into account. 
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1 Introduction 



The theory of Ekpyrotic cosmology 0121121; as well as pre Big Bang scenar- 
ios has emphasized the idea that the Universe did not begin at the Big 
Bang but, rather, had a long prior history. Although the details of this theory 
are far from understood, it is not unreasonable to assume, since the Universe 
had to pass through the Big Bang where densities and temperatures are set 
by the Planck scale, that superstrings play a fundamental role. If this is the 
case, then one expects the natural setting for cosmology to be, not four di- 
mensions, but the higher- dimensional spacetime of string theory. Ekpyrotic 
cosmology introduced the idea that the Big Bang perhaps resulted from the 
catastrophic collision of two brane solitons in this higher dimensional space. 
Colliding branes and the associated bulk space geometry correspond to vac- 
uum solutions of the string theory equations of motion. These vacua may 
solve equations which are valid to some finite order in a string expansion pa- 
rameter, such as five-brane/nine-brane collisions in heterotic M-theory [SlISl 
which are valid to order nf^, or they may be exact conformal field theories, 
which solve the string equations to all orders in the string parameter a'. The 
second type of vacua were emphasized in the Big Crunch/Big Bang j7] real- 
izations of Ekpyrotic theory. Cyclic models jH] are based on Big Crunch/Big 
Bang theories. All of these vacua have in common the property that, prior 
to the Big Bang, the relevant region of spacetime, our past, is contracting 
toward a singular brane collision. In Big Crunch/Big Bang scenarios, the 
Universe then expands outward from this singularity as the Big Bang. The 
singular point is called the "bounce" . 

Of particular importance in all cosmologies, and no less so in Ekpyrotic 
theories, is the origin and momentum spectra of both scalar field and grav- 
itational quantum fluctuations. These are of the utmost importance since 
they produce inhomogeneities in the cosmic microwave background (CMB) 
that have already been observed. These observations are being increasingly 
refined and offer an experimental window into whatever fundamental physics 
is responsible for the present Universe. It was shown in P3|3] that nearly scale 
invariant scalar perturbations are generated in the Big Crunch phase of Ekpy- 
rotic theories prior to the Big Bang. Indeed, a near scale invariant spectrum 
of perturbations is observed in the CMB, but this is in the post Big Bang 
phase. To make contact with these observations, one one must conclusively 
demonstrate that the pre Big Bang fluctuation spectrum in Ekpyrotic theo- 
ries is propagated, nearly unchanged, through the Big Bang. This, despite 
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the fact that Big Crunch/Big Bang geometries are singular at the bounce. 
Although motivated by string theory, almost all previous attempts to ad- 
dress this issue have been carried out within the context of four- dimensional 
toy model geometries, which do not solve the string equations of motion, 
even to lowest order [HI CHI- Recently, this problem was studied in a 
five-dimensional Milne background fTJ ^] , and it was emphasized that the 
five-dimensional structure is important near the Crunch. However, in all 
cases the problem was studied using the techniques of quantum field theory, 
with no string theory input. With no further information at their disposal, 
the authors of these papers had to proceed by matching the relevant pre Big 
Bang and post Big Bang wavefunctions at the bounce, where the geometry 
is highly singular. The choice of boundary conditions is, by itself, conjec- 
tural, being motivated by differing physical arguments. This is made all the 
harder by the singular nature of the geometry at the bounce. The result is 
that some authors claimed that the pre Big Bang scale invariant fluctuations 
are radically altered when they pass through the singularity [H| while others 
claim they are unaltered 112] • The ambiguous nature of these attempts 
was demonstrated in [TU], who showed that, although highly restricted, alter- 
native boundary conditions are possible, leading to the contradictory claims 
in the literature. How, then, can one resolve this ambiguity? 

It seems clear from the previous discussion that string theory should lead 
to a unique solution of this problem. The reason for this is that string theory 
vacua are, in general, globally defined. First of all, their geometric mani- 
folds include both the Big Crunch and Big Bang regions of cosmology, in 
addition to other regions. Secondly, the wavefunctions of these vacua are 
defined everywhere on the geometry. That is, knowing a wavefunction in the 
Big Crunch region, for example, uniquely specifies the wavefunction in the 
Big Bang regime. Clearly, this exactly specifies the boundary conditions for 
the wavefunctions at the singularity, completely resolving the ambiguities 
present in previous work. In this paper, we will show that this is indeed 
the case, at least in the limit of zero string coupling. Working within the 
framework of both supercritical fSlEj and critical type II superstrings^, we 
present a class of Big Crunch/Big Bang cosmological vacua, called "gener- 
alized" Milne orbifolds |15j. The geometry of these vacua includes both the 
past Big Crunch and the future Big Bang regions. There are, in addition, 
four other regions, often called "whiskers": an additional early time region, 

^Note that "type 11" refers to a fermionic string with a chiral GSO projection. 
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an additional late time region and two intermediate regions with closed time- 
like curves. The intermediate regions connect to the Big Crunch/Big Bang 
regions at the bounce singularity. These vacua are exact superconformal 
field theories and our results are valid to all orders in the string worldsheet 
parameter a' ^3]. In this paper, we will present our discussion within the 
context of supercritical string theory, since the corresponding spacetimes 
are manifestly homogeneous and isotropic. However, we can show that the 
results are identical for certain solutions of critical string theory. These ten- 
dimensional vacua contain, in addition to the generalized Milne directions, 
two non-compact spatial directions described by the two-dimensional "cigar" 
conformal field theory ^B]; which we take to be two of the three large space 
directions. These vacua are not homogeneous and isotropic, but approach 
these properties, as closely as one likes, if a parameter is taken to be small. 
The Milne orbifold can be obtained as a specific limit of these generalized 
Milne orbifolds. Working within this context, we find that the wavefunctions 
are, indeed, globally defined. This follows from the fact that the generalized 
Milne vacua all descend from a "covering space" by the process of cosetting 
out a gauge action and orbifolding. The invariant globally defined wavefunc- 
tions on the covering space then descend to globally defined wavefunctions on 
the generalized Milne orbifolds. Therefore, to all orders in a', the boundary 
condition ambiguities inherent in previous work have been resolved. 

What, then, are the results for the fluctuation spectrum? We find that 
the fluctuation spectrum in the Big Crunch regime is, in general, altered by 
its passage through the Big Crunch/Big Bang singularity. The change in 
the spectrum can be calculated unambiguously at any time after the bounce. 
In the far future, it can be expressed by an explicit momentum and time- 
dependent function A{k,t), which multiplies the early time pre Big Bang 
fluctuation spectrum. In the Milne orbifold limit, we flnd that A — >■ 1 and, 
hence, the Big Crunch fluctuation spectrum is preserved as it passes through 
the singularity. This proves the conjecture flrst introduced in [3] and shown 
within a flve-dimensional field theory context in [12|. If this result survives 
corrections due to gravitational backreaction and stringy effects, which we 
will comment on at the end of this Introduction, it means that the fluctuation 
spectrum of Ekpyrotic cosmology may well be consistent with observation. 
However, for generalized Milne orbifolds A is not unity and the fluctua- 
tion spectrum changes as it passes through the bounce singularity. We will 
show that this change is entirely due to the existence of the whisker regions. 
Speciflcally, the quantum mechanics of these stringy regions is inextricably 
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linked to the quantum mechanics of the Big Crunch/Big Bang geometry. 
The point of hnkagc is at the bounce singularity, where the whiskers and the 
Big Crunch/Big Bang regions touch. The result is that there is, in general, 
particle production in the future Big Bang region, even though the vacuum 
in the Big Crunch was empty. This particle creation affects all of the correla- 
tion functions in the future. In particular, it affects the two-point correlation 
function from which A is calculated. We find that A is an explicit function 
of momentum, with both time-independent and time-dependent components. 
This calculation exposes a subtlety that arises whenever the geometry has 
whiskers of the type discussed in this paper. That is, that the early time in- 
vacuum is ambiguous. If we define this vacuum to be such that an observer 
in the Big Crunch regime sees no particles, then we find the associated in- 
vacuum in the early time whisker is not completely fixed. This ambiguity in 
the in-vacuum can be parameterized by two constrained complex functions 
j{k) and j{k). The parameter j{k) explicitly enters the expression for A. 

— * — * 

One natural choice turns out to correspond to ^{k) — l,7(/c) = and leads 
to particle creation and a deviation of A from unity. However, as 7 0, 
the particle production in the Big Bang region decreases to zero and A ^ 1. 
Therefore, for a large choice of in- vacua, the change in the fluctuation spec- 
trum by the singularity is small. This might open the possibility of measuring 
string effects, namely the existence of whisker regions connected to the Big 
Crunch/Big Bang geometry, as small momentum and time-dependent devia- 
tions from the scale invariance of the inhomogeneities in the CMB. We hope 
to discuss this further elsewhere. 

The hnkage of the quantum mechanics of the Big Crunch/Big Bang and 
whisker regions has a second, unanticipated effect. As viewed in the Hilbert 
space of the complete quantum mechanics, the in-vacuum is a pure state. 
However, we find that it is "entangled" , that is, contains correlations between 
the late time Big Bang and whisker regions. Therefore, tracing its density 
matrix over the states of the unobserved whisker, one obtains a non-trivial 
density matrix in the observable Big Bang region. That is, an observer in 
the Big Bang region, with no access to information about the whiskers, finds 
himself in a mixed state. This state has "entanglement" entropy, which man- 
ifests itself explicitly in the expression for A. If the entropy were zero, then 
the expression for A should be compatible with a Bogolubov transformation 
linking the "in" and "out" states of the Big Crunch and Big Bang regions. 
However, non-vanishing entanglement entropy will ruin this compatibility. 
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We find that our explicit expression for A indicates non-vanisliing entangle- 
ment entropy, except in the limit 7 — where the whisker effects decouple. 
It follows that, not only are the changes in the spectrum calculable, but they 
explicitly exhibit entropy induced by the existence of the stringy whisker re- 
gions. This is not dissimilar to recent discussions of entangled states within 
the context of BTZ black holes ^7], see also |18j . 

We have presented the results of this paper mostly in the context of 
Ekpyrotic Big Crunch/Big Bang transitions. Indeed, Ekpyrotic theory has 
inspired much of the recent interest in Big Crunch/Big Bang singularities in 
string theory [71 and presents a context where it is particularly important 
to know what happens to fluctuations at such a singularity. However, Big 
Crunch/Big Bang transitions are of more general interest and are especially 
appealing in string theory, where observables are usually defined using simple 
asymptotic regions. See, for instance, Therefore, we believe that the 
results of this paper, such as information loss due to the existence of whisker 
regions, might be of interest in more general cosmological scenarios. 

Specifically, in this paper we do the following. Section 2 is devoted to 
a discussion of a specific four-dimensional quantum field theory involving 
gravity and three scalar fields. In subsection 2.1, we present a cosmological 
background solution of this theory and explore various physical properties. 
It is shown that there are two independent branches to this solution, one 
representing a Big Crunch universe, evolving from the past and ending in 
a singularity, and the other a Big Bang universe, beginning in a singularity 
and evolving into the future. Scalar field quantum fluctuations, in each of 
the Big Crunch and Big Bang regions, are discussed in subsection 2.2. Both 
the in-vacuum and the out-vacuum are defined, and we explicitly compute 
the scalar two-point correlation function with respect to each of these vacua. 
In subsection 2.3, we relate both the classical and quantum theories of the 
Big Crunch and Big Bang regions by connecting them at the singularity. 
Within this context, we compute the two-point correlation function with 
respect to the in-vacuum in both regions and compare them. We find that 
the scalar fluctuation spectrum is potentially altered when it passes from the 
past Big Crunch region through the bounce singularity. In the far future, 
the change in the spectrum can be expressed by a momentum and time- 
dependent function. We compute this function explicitly and show that it 
depends on the Bogolubov coefficients relating the "in" and the "out" states 
of the Hilbert space. The meaning of these coefficients, and how one should 
compute them, is discussed. Finally, in subsection 2.4, we show that our 
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four-dimensional quantum field theory is the low energy effective theory for 
supercritical type II strings. 

Section 3 is devoted to superstring cosmology within the context of type 
II supercritical string theories. In subsection 3.1, we review the "general- 
ized" Milne orbifold solutions of these theories first presented in ^]. These 
solutions are exact superconformal field theories. The various regions of the 
associated spacetimes are discussed in detail, including the additional whisker 
regions. It is shown that, at low energy, these string vacua give rise to the 
Big Crunch/Big Bang theories introduced in section 2. We study the scalar 
fluctuation spectrum in subsection 3.2. The structure of generalized Milne 
orbifolds as the coset and orbifold of a PSL{2, IR) covering space is reviewed 
and the relationship of these vacua to the Milne orbifold is discussed. In sub- 
section 3.2.1, a basis of wavefunctions, each defined on every region of the 
spacetime, is presented and studied. We quantize the scalar fluctuations in 
subsection 3.2.2. This is accomplished by expanding the scalar fluctuations 
in this basis and canonically quantizing the coefficients. We define two nat- 
ural vacua, the in-vacuum and the out- vacuum, and compute the Bogolubov 
coefficients relating them ^3]. Using these results, we compute the particle 
production in the future regions and the scalar two-point correlation function 
with respect to the in-vacuum. This result is valid in any region, allowing us 
to compare the spectrum in the future Big Bang region with the early time 
spectrum during the Big Crunch. We find that, generically, the spectrum 
is altered. In the far future, we can express the change in the spectrum in 
terms of a function A, which is momentum and time- dependent. We com- 
pute this function explicitly. In subsection 3.2.3, we show that there is, in 
fact, a family of in- vacua, which we specify with two constrained complex 
functions •y^k) and 7(fc). We repeat the calculation of the Bogolubov coef- 
ficients, particle production, scalar two-point correlation function and A in 
this context. Again, in general, the fluctuation spectrum is changed as it 
passes through the bounce singularity. However, in the limit that 7 — 0, 
A approaches unity and the spectrum is conserved. In this limit, A can be 
expressed in terms of the pure Big Crunch/Big Bang Bogolubov coefficients. 
However, for any finite 7, this is no longer the case. This is explained in sub- 
section 3.3, where it is shown that the generalized in-vacuum is an entangled 
state in the future. This implies that, from the point of view of an observer in 
the Big Bang region, this vacuum has entanglement entropy. This is equiv- 
alent to quantum mechanical "information loss" into the whisker regions. It 
is this information loss that obstructs writing A in terms of the pure Big 
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Crunch/Big Bang Bogolubov coefficients. In subsection 3.4, we discuss the 
limit of our superconformal vacua to the Milne orbifold. It is shown that, 
in this limit, the factor A goes to unity for any 7 in- vacuum. This proves 
that, in the Milne orbifold cosmology, the fluctuation spectrum is unaltered 
by the bounce singularity. In subsection 3.5, we include some comments on 
backreaction. Finally, in Appendices A, B, and C, we outline the theory 
of global wavefunctions, present the expressions for a basis of these wave- 
functions in all regions of the spacetime and discuss the Milne limit of these 
wavefunctions respectively. 

Before proceeding, we would like to make several important comments. 
First, note that the generalized Milne orbifolds have a smooth circle in the 
spatial Milne-direction. However, it is not hard to show that all of the results 
of this paper will be unchanged if we allow a further IL2 orbifolding in this 
direction. In this case, this smooth circle becomes a finite interval, with new 
twisted sector states appearing on each of the two boundaries. This new Z2 
orbifolded vacuum is again an exact superconformal field theory, which is 
closer in spirit to the notion of colliding branes. 

As we have stated previously, our results are computed with respect to 
exact superconformal field theories, so a' corrections are under control. How- 
ever, it is well-known that these theories suffer from severe backreaction prob- 
lems, which could conceivably modify our results. Indeed, string perturbation 
theory gives rise not only to an expansion in a' (higher derivative corrections 
to the classical action), but also to an expansion in the string coupling Qs 
(quantum corrections). The results of ^1120] on global wavefunctions refer 
only to the first term in the gs expansion. It has, in fact, been shown in 
following [22], that classical string scattering amphtudes exhibit divergences 
associated with a large backreaction of the spacetime geometry to small per- 
turbations. For a non-perturbative manifestation of gravitational instability, 
see Large gravitational backreaction (or the absence thereof) in this and 
related models j20^i211l2S] was also studied in [22] • The tree- level divergences 
of [211 122] indicate a breakdown of string perturbation theory. That is, it is 
not consistent to ignore gg corrections if gs is small but non-zero. For a non- 
technical discussion of this, see [27]. Therefore, there is at present no fully 
controlled computational framework, and it is conceivable that our detailed 
results will receive significant corrections when backreaction is taken into ac- 
count. One may hope that at least certain important qualitative features, 
such as the peculiar causal structure of the string solutions and its effect on 
four-dimensional cosmology, will survive gs corrections. However, it will take 
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significant advances in string theory to either estabhsh or refute this. 

It has been suggested, see for instance j2Hl, that string winding modes 
that become hght near the bounce might play a role in resolving the singu- 
larity in string theory. In the case of the Milne orbifold, the description of 
these winding modes is somewhat complicated, because the size of the Milne 
circle grows without bound away from the singularity. See [29j for a very 
recent discussion and ^D] for earlier work. However, in the generalized Milne 
orbifold the radius of the Milne circle approaches a constant asymptotically, 
and the vertex operators for the winding modes are explicitly known ^ISJ. In 
the present paper, we will ignore winding modes, except to note that they 
become light near the Big Crunch/Big Bang singularity and invalidate a four- 
dimensional description there. They are heavy away from the bounce and, 
thus, do not appear in the effective action describing our four-dimensional 
cosmology. However, we should stress that the computations we do for the 
generalized Milne orbifold can easily be extended to include winding modes 
|15j . Of course, the most interesting effects of winding modes should involve 
string interactions, and these have not yet been computed for the generalized 
Milne orbifold. 

In string theory, the observables are S-matrix elements. At zero string 
coupling, the regime in which we will be working, they are determined by a 
Bogolubov matrix. Our strategy is to compute this matrix using the glob- 
ally defined vertex operators jT^l 1201 • More specifically, we calculate the 
Bogolubov matrix to leading order in a' from the associated globally defined 
wavef unctions. However, as was mentioned in |[15j, the higher a' corrections 
modify this matrix in a very trivial way, simply multiplying some of the 
entries by a phase. Therefore, the particle creation rates we compute are 
exact to all orders in a'. In this paper, we will go a step beyond computing 
S-matrix elements and compute correlation functions at a fixed finite time, 
using the global wavefunctions we obtained from string theory. 

For additional recent string theory work related to Big Crunch/Big Bang 
singularities, see pij . 

2 Four- Dimensional Cosmology 

In this section, we will explore the cosmological properties of a specific four- 
dimensional quantum field theory coupled to gravity. This theory consists 
of three scalar fields, denoted by ctt, cr^ and respectively, coupled to the 
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usual Einstein gravity. The action for this theory is given by 




where 

C^R- \g^''d,<t>d,(t> - \g^''d,aTd,aT - \g^''d, andean - AQ'et (2) 

There are two dimensionful parameters in this action, Newton's constant 
^4/87? (^4 having the dimension of a length) and a mass scale Q, which we 
will think of as small compared to the momenta of interest. 



2.1 The Background Spacetime 

We will be interested in cosmological solutions of the equations of motion 
of this theory that are spatially homogeneous and isotropic, that is, of the 
Friedman-Robinson- Walker (FRW) type. We find the following general class 
of such solutions. The scalar fields are independent of all spatial coordinates 
and evolve as follows. The field (Tt simply vanishes, 

ot = 0, (3) 

whereas 

ctr = log I tanh(Qt) I (4) 

and 

= 20o-log|sinh(2Qt)|. (5) 

Here 0o is an arbitrary integration constant. These fields drive a time- varying 
metric which is independent of all spatial coordinates and given by 

g^iv = aitfri^^, (6) 

where 

^2 ^g-2<^o|sinh(2gt)| (7) 

and rjni, is the metric of flat Minkowski space. It follows from the form of the 
metric that t is conformal time. 

What are the cosmological properties of this class of solutions? Perhaps 
the most salient feature is that the conformal factor of the metric, a(i)^, is 
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monotonically decreasing in the time interval — oo < t < and monotonically 
increasing for < t < +00. We will refer to the spacetime in the negative 
time interval as Region I and the spacetime for the positive time interval as 
Region II. Let us be more specific about the geometry of these regions. To 
make contact with the conventional analysis of FRW cosmologies, we will 
first compute the equation of state 

w = ^, (8) 

where P is the pressure and p the energy density. These quantities are defined 
by the Einstein equations 

= '^a^p (9) 



and 

W = - 

6 

where 



U' = -'^a\p + W), (10) 



and the prime denotes differentiation with respect to conformal time t. Us- 
ing dZj), we find from Q and ^ that 



^ kI |sinh3(2gt)| ^ ' 



and 

p ^ Q^^2<^„ / 4-cosh^(2gt) \ 
k| V |sinh^(2gt)| ) 

respectively. It follows that 



(13) 



3 V cosh2(2Qt) I ^ ' 

The first thing to notice is that w is not constant in time. In the far past 
and future ^ 

w — »• — , t — »• ±cxo. (15) 

3 

As the cosmology evolves, one finds that 

w = 0, to = ±7^ log(2 + VS) ^ cosh(2gto) = 2 (16) 
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and 

w = -, ti/3 = ±^ log(l + V2)^\ sinh(2gti/3) I = 1. (17) 

Finally, as one approaches the origin, either from negative time or positive 
time, 

t-^±0. (18) 

It follows that 

- ^ 1. \t\ « g. (19) 

What is the interpretation of this somewhat unusual FRW cosmology? To 
elucidate this, note that in a scalar dominated phase the energy density and 
pressure are given by 

where ipi, i = 1, 2, 3 are the fields cxr, cfr and (j) respectively and V is the total 
potential energy. Using (IS)), © and it is easy to show that p and 
P in ()2U|) are identical to expressions ()12|) and ()13|). It follows that in both 
Regions I and II our cosmology is scalar dominated with the equation of state 
w changing as the scalar fields evolve. In the conformal time coordinate t, 
only in the regions described by (fTH|) and (fTIHl does w behave approximately 
as a constant, w = —1/3 and w = 1 respectively. Note that, in the latter 
time region, the conformal factor a{t) in becomes 

9 = ^, (21) 



'e2<^o/2Q' ' ^ 2 

which is consistent with constant w = 1. 

A second important feature of FRW cosmologies is the Hubble parameter 

H = -, (22) 

and its inverse, Rh = \H\~^, the Hubble radius. Using ((Tj), one finds that 

H = Qef^ coM^Qt) .23) 

sinh(2Qt)|sinh(2Qt)|^ 
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Again, note that H and, hence, Rh are not constants but, rather, evolve 
with conformal time t. Let us focus on the behaviour of the Hubble radius. 
In the far past and future 

Rh — > oo, t — ^ ±00. (24) 

The Hubble radius then changes monotonically, taking the value, for example, 

6"''^° 1 r- 

Rh = ^=-, ti/3 = ±777^ log(l + V2) (25) 



before behaving as 



' {I -QH' + ...), t-^±0 (26) 

as the time approaches the origin. Hence, 

Rn - \t\ « (27) 

which vanishes at t = 0. Note that this behaviour is consistent with the 
monotonic evolution of the conformal factor and is valid in both Region I 
and II. From ((7j) and (j27j) it is clear that the Hubble radius goes to zero more 
quickly than the physical wavelength of an excitation, which is proportional 
to a. Therefore, all modes are "frozen" outside of the horizon at the Big 
Crunch/Big Bang transition. Similarly, we conclude from ((7j) and (j^^ that 
at very early and very late times, modes are inside or outside the horizon 
depending on whether their comoving momentum k satisfies k"^ > or 
fc^ < Q^, respectively. In the following, we will focus on fluctuations with 
k^ ^ Q^. Thus, these modes start out inside the horizon at very early times 
and freeze as they approach the Big Crunch. Similarly, in the Big Bang 
region of spacetime they start out frozen near the Big Bang and enter the 
horizon at some later time. 

A third important quantity to consider is the scalar curvature, R. In 
conformal time, the scalar curvature in FRW spacetimes is given by 

(28) 

It then follows from ((7j) that 

sinh^(2Qt) - 1 

« = ^«'*WW' '''' 
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an expression that is valid in both Regions I and II. Key values of the scalar 
curvature occur at precisely the same times, namely ±00, to, ti/3 and ±0 
defined in fll6|) and ()17|) respectively, that indicated the behaviour of the 
equation of state w. In the far past and future 

R — ^ 0, t — > ±00. (30) 

As \t\ decreases from infinity toward zero, R has the following properties. 
To begin with, R is positive and increasing until l^ol = 2^ log(2 + after 
which R monotonically decreases. The scalar curvature vanishes at 

i? = 0, ti/3 = ±^log(l + v^). (31) 
For smaller values of R is negative. As one approaches the origin, 

i?— > -7e'*°^(l-6QV + ■■■), t-^±0 (32) 



which diverges as 



4 t 



R = -^e^^o^, \t\ < (33) 



The above results allow us to give a concise description of our cosmological 
solution. First consider Region I, corresponding to the negative time interval 
— 00 < t < 0. The associated geometry is that of a spatially homogeneous and 
isotropic FRW spacetime. In the distant past, the manifold has vanishing 
scalar curvature and a divergent Hubble radius. As time progresses, the 
scalar curvature first grows positively, reaches a maximum and then begins 
to decrease, vanishing at a finite time ^1/3. Henceforth, the curvature is 
negative, diverging as as t approaches the origin t = 0. Throughout 
the entire time interval —00 < t < 0, the Hubble radius is monotonically 
shrinking from infinity to zero. Both the vanishing of the Hubble radius 
and, particularly, the divergence of the scalar curvature as t — > 0, tells us 
that Region I terminates abruptly at t = 0. Region I, therefore, is a classic 
example of what is called a "Big Crunch" cosmology. Region II, on the 
other hand, corresponding to the positive time interval < t < +00, is the 
exact mirror image of Region I in the time direction. That is. Region II is a 
spatially homogeneous and isotropic FRW spacetime that starts abruptly at 
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t = with vanishing Hubble radius and negatively infinite scalar curvature 
and then expands outward. The scalar curvature increases from negative to 
zero to positive, reaches a maximum and then decreases to zero as t — > +00. 
During the entire time interval < t < +00, the Hubble radius monotonically 
increases from zero to infinity. Region II, therefore, is a classic example of 
a "Big Bang" cosmology. It is essential to note that in general relativity, 
because of the curvature singularity at t = 0, there is no relationship between 
Region I and Region II, each representing an independent cosmology. 

2.2 Fluctuations 

We now turn to a discussion of quantum fiuctuations of the scalar fields on 
the Big Bang/Big Crunch geometries presented above. To do this, we must 
first expand (Tt, ctr and (j) around their classical values, which we now denote 
as (o"t), {cr) and (0), given in (Hj) and (jS)) respectively. That is, 

arp = {(Xt) + S(Tt, an = (a/j,) + 6(Tr, = (0) + Sc/). (34) 

Inserting each of these fields into its equation of motion, using expressions ^ 
and ((Tj) for the background metric and assuming the fiuctuations are of the 
form 

= 5T(t)e^^^ 6aR = 6R{t)e'''-^ , 50 = 5$(t)e^^^ (35) 
we find that the fiuctuation 6T is a solution of 

6T" + 2Q coth{2Qt)ST' + P6T = 0, (36) 
6R solves the same equation 

6R" + 2Q coth{2Qt)SR' + PSR = (37) 

whereas 

5$" + 2Q coth{2Qt)6^' + {P + AQ'^)5^, 5$ < 1. (38) 

Note that the 5$ fiuctuations will satisfy the same equation as 5T and 5R 
for momenta 

P > 4g2. (39) 
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Henceforth, we will restrict our discussion to this momentum regime. Since, 
in this case, all three fluctuations are specified by the same equation, we will 
simply focus on one of them, which we choose to be 5T. 

Let us search for solutions of the 5T fiuctuation equation (jHBj) . To do 
this, we must first specify the region of spacetime in which we want to work. 
Begin by considering Region I, with negative conformal time in the interval 
— oo < t < 0. For very early times, (j36|) simplifies to 

5T" - 2Q6T' + P6T = 0. (40) 

It is easy to see that this has plane wave solutions of the form 

ST^ = Qe^^e^^^s*, (41) 

where 



E^^ = ^k^ - Q2 (42) 

is the energy associated with momentum k and is a normalization con- 
stant. Note that is a positive real number in the momentum regime ()39p 
in which we are working. The normalization constant can be determined 
using the scalar product 



/ 0i9m0;v^c/S^ (43) 



/s 

where E is a space-like three-surface, is the volume element on that 
surface, 

dY.^ = 5\^J^dx (44) 

and 0j, i = 1,2 are any two scalar functions. Using the metric given in ^ 
and ©, expression ()41|) for ST and ()43p . we find that 

(5T^e'^-", ST^e'^'-"^) = ±(27r)=^5(A^ - k')\C^\\-^^° E,;. (45) 

Note that the Klein-Gordon norm obtained from ()43|1 can be either positive or 
negative depending on the frequency of the scalar function. Setting the right 
hand side of ()45|) equal to ±{27r)^S{k — k'), it follows that the normalization 
constant, up to a phase, is given by 

p4'o 

Q = 4^- (46) 
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Therefore, the normahzed asymptotic plane wave solutions of (|4Uj) are given 
by 



ST^ = -^e^'e^'^x\ t -oo. (47) 



One can, in fact, solve the 6T fluctuation equation (j36p for any value of t in 
Region I. The result is found to be 

Ajp<l>o 1 

= -j={-zyF{-j, -r, -2r, -), (48) 

where 

J = -^ + ^^, ^ = -sinh2(gt) (49) 

and F{a, b; c; x) is the hypergeometric function 2F1 (see Appendix B for its 
definition and some properties). In addition, the hermitian conjugate 

5Tr = 6Tf* (50) 

is an independent solution of (jHUjl . Using the facts that 

z — > — , F{-j, -j; -2j; -) — ^ 1 as t — > -00, (51) 



we see that fl48|l and ()5()|1 approach the plane waves solutions ()47j) in the far 
past. One can show that (j48|) and (j5(J|) diverge logarithmically at t = 0. 

Combining these results with the first expression in (j35|) . we see that any 
fluctuation Sax in Region I can be written as 



^^(a/5T^ (t)e^^-^ + a/*6T^ {t)e-^''-^), (52) 

where a^^ are arbitrary complex coefficients. Note that the first function in 
this expansion, 5T^{t)e^^'^ , corresponds to a pure positive frequency plane 

wave in the far past. Similarly, in this limit, 5T^ {t)e~'^^'^ becomes a negative 
frequency plane wave. 

Thus far, our discussion of the fiuctuations 5aT^ has been strictly classical. 
However, the theory can be easily quantized by demanding that 5(Jt^ and, 
hence, the coefficients a^^ be operators in a Hilbert space. For simplicity. 
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we will continue to denote these operators by 5aT^ and a^^, suppressing the 
usual "hat" notation. The quantization will be canonical if we assume that 

%-'^'] = - k% [a/, a,-/] = [agt/^ ^ .f/] = q. (53) 

The vacuum state of the quantum theory is then defined as the normalized 
state |0)i„ satisfying 

a/|0)i„ = (54) 

for all momenta k. There are many objects that can now be discussed in this 
context. In this paper, we will focus primarily on the two-point correlation 
function 

i„(0|(5a/(t, x)6aT\t, x + r)|0)i„. (55) 

Using (fS^. and we find that this function is independent of x and 
given by 



1^1 1^ 



27r 

In the following, we will always set 7^= and consider 

,„(0|5a/(t,x)».„ = j (0|5Vr, (57) 

since this is sufficient for determining the fluctuation spectrum. The result for 
general r can be obtained by multiplying the integrand by sin(|/c||?^)/|/c||?^. 

Equation (|57p can be computed for any time t in Region I using the 
expression given in (j48|) . However, for our purposes, it is most illuminating 
to evaluate it in the distant past. Inserting expression ()47|) . we find that the 
correlation function (|37|) becomes 

,„(0|5a/(t,f)2|0),„ = m J (58) 

where 

f{t) = 2e2<^°e2«*. (59) 
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For ^ Q^, the momentum regime (j39|) in which we are working, expres- 
sion becomes, to next to leading order, 

.„(0|fa.'«, 2fm.„ = /(i) / l^j^ (1 + II?) . (60) 

The momentum dependence of the first term of the integrand is simply that 
of zero-point fluctuations in Minkowski space. On the other hand, the sec- 
ond term corresponds precisely to a scale invariant fluctuation spectrum. 
Note, however, that since ^ Q^, the second term is subdominant to the 
Minkowski fluctuations. The time- dependent factor f{t) is not canonical. 
To understand its origin, we note that the kinetic energy term for cxr in the 
Lagrangian Q is not canonically normalized in the gravitational background 
given in © and (|7j). This kinetic energy term can be canonically normalized 
by defining a new scalar field as 

= e-<^« I sinh(2Qt) | ^/V^. (61) 

Note that in the far past this expression becomes 

St' = /(t)-'/VT, t -oo, (62) 

where f{t) is given in It follows from this and (jUUI) that at early times 

, ,c- ^ f d^k I (l / 2\ , , 

,„(0|ffi/(«.f)^|0>,„ = / (5 + ^] ■ (63) 

We conclude that the factor /(t) in the correlation function (jHn|) is simply a 
scale factor that can be absorbed by a field redefinition or not, depending on 
taste. This concludes our analysis of quantum fluctuations in Region I. 

We now want to discuss the quantum fluctuations of boT in Region II, 
that is, for conformal time in the positive interval < t < +00. Note that the 
fluctuation equation (j36j) is independent of which region is being considered. 
It follows that the analysis of quantum fluctuations in Region II is essentially 
identical to that in Region I. For this reason, we will simply present our 
results. To begin with, the fluctuations bTi and (5Tr given in PHj) and (jSOJ 
respectively remain a complete set of solutions of equation (p?H|) . It follows 
that any quantum fluctuation in Region II can be written as 

^<yT' = J jS^M^'^n^^y'' + <^"'^T^iiy~''% (64) 
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where and a,% satisfy the canonical commutation relations 

[a/^ a,-/^t] = - k% [a/^ ag'^] = [a/^\ a,-/^t] = q. (65) 

The vacuum state is then defined as the normalized state \0)out satisfying 

a/'\0)out = (66) 

for all momenta k. Again, there are many objects that one may wish to 
compute at this point. For example, in analogy with Region I, we find that 
in the distant future 

ouMScTT'^t, xflOU = 9{t) / + ^) , (67) 

where 

git) = 2e^'^«e-2^*. (68) 

As discussed previously, the factor of g(t) arises from the non-canonical nor- 
malization of the (Jt kinetic energy term in ((21) with respect to the geometric 
background © and ((7|). Proper normalization of this term can be achieved 
by defining 

j^^n = g-<^o I sinh(2Qt) | Vt'' (69) 

for any positive conformal time t. Note that in the far future this relation 
becomes 

E/^ = gitr'/'ar'', t — +oo, (70) 

where g{t) is given in ()68p. In terms of this canonically normalized scalar 
field, the fiuctuation spectrum ^7\i becomes 




Therefore, the g{t) factor in the two-point correlation function ()67j] is simply 
a scale factor. It can be absorbed or not, depending on taste. This concludes 
our analysis of quantum fiuctuations in Region II. 

As stated earlier, because of the curvature singularity at t = 0, there is 
no classical relationship between Region I and Region II. Each represents an 
independent cosmology. The same is true for the quantum fiuctuations that 
we have just discussed. The creation operators a^^^ acting on the vacuum 
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|0)j„ create a Hilbert space of states Ti^ representing the quantum theory in 
the Big Crunch geometry of Region 1. Similarly, the operators a^^^^ acting on 
|0)out create a Hilbert space 7i^^ representing the quantum theory of the Big 
Bang geometry of Region II. A priori, there is absolutely no relation between 
and n^^. 



2.3 Relating the Big Bang to the Big Crunch 

It is the thesis of Ekpyrotic cosmology that the Universe did not begin at the 
Big Bang. Rather, it had a prior history, connected to our present geometry 
via some catastrophic event. In the Big Crunch/Big Bang versions of Ekpy- 
rotic theory, this catastrophic event is a spacetime singularity at t = 0. This 
singularity is of the type found in the curvature scalar in Region I and Region 
II as if: —i> —0 and t — >■ +0 respectively. We will, therefore, within the context 
of the theory described by (|21), construct a Big Crunch/Big Bang scenario 
by connecting Region I and Region II classically at the singular point t = 0. 
Having done this, one must also specify a relation between the quantum the- 
ories on these two regions. The most naive approach, and the one we will 
adopt in this section, is to identify the two Hilbert spaces. That is, assume 
that 

= = n. (72) 

One consequence of this is that the creation/annihilation operators a^:^ ,a^^'^ 
and a^^^a^:^^'^ all act on the same Hilbert space Ti. In "normal" quantum 
field theory, that is, when there is no geometric singularity or event horizon 
separating the past from the future, the "out" creation/annihilation oper- 
ators are linearly related to the "in" creation/annihilation operators via a 
so-called Bogolubov transformation. We will assume that the same is true 
in our theory, despite the existence of a curvature singularity at t = 0. That 
is, we postulate that 




X*{k) Y*{k) 
Y{~k) X{-k) 




(73) 



The complex Bogolubov coefficients X and Y are not completely indepen- 
dent. They are constrained by the requirement that the Region I and Region 
II creation/annihilation operators continue to satisfy the canonical commu- 
tation relations (j53|) and (j65j) respectively. It follows that 

|Xp-|F|2 = l. (74) 
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This is the only constraint on these coefficients. However, for simphcity, we 
will further assume that 

X{k)=X{-k), Y{k)=Y{-k). (75) 

Relaxing these assumptions will not change any of our conclusions. 

Having postulated the relations (f?^ and (fTSj) . one can now compute 
correlation functions that were not defined separately in Region I and Region 
11. Specifically, we want to calculate the two-point function 

,„(0|5a/'(t,x)2|0),„, (76) 

where Sax^^ is the Region II field operator given in (jMjl . whereas |0)j„ 
is the Region I vacuum defined in (j54j) . This is easily accomplished us- 
ing dnni), (El, (El, dZSl) and dZl. The result is 

in{0\5aT'\t,xf\0)in = 

= /(0j((l + 2|rp)|5T+p + XWf + X*y*5Tr2), (77) 

with the fluctuations 6Tf and 6T^ given by (jIH|) and (EUj) respectively. Of 
particular physical importance is the form of this correlation function in the 
distant future. As t ^ +oo, expression (f77j) becomes 



in 




where the scale factor g(t) is given in and 

A{k, t) = l + 2|F|^ + XYe-^'^>^' + X*Y*e^'^^'. (79) 

Written in terms of the scalar field 'Zt^' defined in (EHl), and using (fTOj). this 
simplifies to 

..(0|5S/^(t,x)^|0).„ = / (1 + ^) Aikt). (80) 

We want to compare these results to another correlation function, namely 

in{0\SaT\t,xf\0)in. (81) 
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Note that, in this case, both the field operator 5aT^ and the vacuum |0)j„ are 
the Region I quantities defined in (|52p and (|54j) respectively. For arbitrary 
time t, this two-point function was evaluated in ()56|) and its functional form 
in the limit t — *• —oo presented in fl60|) . Finally, in terms of the scalar field 
defined in (jHH), and using the t ^ —oo limit of this correlation 
function becomes 

.„WE/(M^^|0>.„ = /|^5^(i + |g). (82) 

Although these expressions were first derived strictly within the context of 
Region I, they remain valid for the complete Hilbert space Ti. First note 
that, for arbitrary conformal time t, correlation functions (j77|) and (j56p are 
identical only if the complex Bogolubov coefficient F = 0. This remains 
true when comparing the t — >• +00 limit of m(0|5(TT'^^(t, xY\^)in given in (fTSjl 
to the t —00 limit of i„(0|(5(jT^(t, x)^|0)j„ presented in (jHH). The simplest 
comparison can be made between the future correlation function (j80|l and the 
past correlation function ()82|1 . since the irrelevant scale factors have been 
removed. As stated earlier, the form of the argument of the momentum 
integral in (jH^ is that of Minkowski fiuctuations plus a subdominant scale 
invariant contribution in the far past. However, the A{k,t) factor in ()80|] 
potentially modifies the fiuctuation spectrum in the far future. If coefficient 
Y = 0, then it follows from (j79p that A = 1 and one again obtains the 
spectrum (I82|) . However, if F 7^ 0, then A 7^ 1 and the fiuctuation spectrum 
is modified. What are the physical consequences of this? 

A central assertion of Ekpyrotic Big Crunch/Big Bang theories is that a 
scale invariant fiuctuation spectrum is generated in the Big Crunch geome- 
try prior to the singularity which is then transmitted, without modification, 
to the Big Bang geometry after the singularity. It is these scale invariant 
fiuctuations that are assumed to account for the observed fiuctuations in 
the cosmic microwave background. But is this true? Or is the fiuctuation 
spectrum modified by the presence of the singularity? There has been con- 
siderable controversy regarding this, with some authors concluding that the 
spectrum is transmitted unchanged, some authors claiming it is greatly mod- 
ified and further literature showing that this question is ambiguous as posed, 
requiring more physics input to uniquely resolve it. All of this literature has 
attempted to confront this issue by imposing explicit boundary conditions to 
match the incoming and outgoing wavefunctions at, or near, the singularity. 
However, attempting to study the vicinity of a singularity is difficult since 
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one expects short distance effects to greatly modify the geometry and the 
physics. This is the source of the ambiguities. Our point of view is differ- 
ent. We see that the question of the persistence of the spectrum through 
the singularity is precisely expressed by whether or not the Bogolubov coef- 
ficient Y vanishes. If it vanishes, the spectrum is transmitted through the 
singularity into the future. If it does not vanish, the spectrum is modified 
after the Big Bang. We conclude that to study this problem, one must, un- 
ambiguously, compute the Bogolubov coefficient Y. But how? One could, of 
course, attempt to compute Y by matching the wavefunctions at the singu- 
larity. But, as we have said, this process would be ambiguous. A far more 
concrete way would be to compute Y directly from string theory, where at 
least for certain types of singularities global wavefunctions can be unambigu- 
ously defined. Given such globally defined wavefunctions, one can perform 
calculations in the asymptotic past and future, far away from the singularity 
at t = 0. This is the approach we will follow in the remainder of this paper. 
However, to compute the Bogolubov coefficient Y in this manner, it is nec- 
essary to demonstrate that the field theory defined by © is, in fact, the low 
energy four- dimensional effective theory of some specific string theory. This 
is indeed the case, as we will now show. 

2.4 Four-Dimensional Cosmology from String Theory 

In ((i+l)-dimensional spacetime, string theory gives rise to a classical effective 
action of the form 

Sd+i = I f^'+'a;v^e-2*(i?,+i + Ag^'dmi^ - 2A), (83) 

where /, J = 0, 1, . . . , (i, gu and $ are the [d + l)-dimensional string frame 
metric and dilaton respectively and all other massless and massive string 
modes have been set to zero. A positive tree level cosmological constant 
A > will arise, for example, in supercritical type II string theories in 

d+ 1 = 26,42,58,... (84) 

dimensions |14j . 

Let us now assume the existence of solutions of the string equations of 
motion with metrics of the form 

gijdx^dx^ = g' dx^dx" + e'^'^Tdabdx^dx^ + e^'^^^dx'^f , (85) 
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where /i, z/ = 0, 1, 2, 3, indices a,b = 4, . . . ,d — 1 and g'^^^, a'rp and cr^ are 
all functions of the four- dimensional coordinates x^, /i = 0, 1,2,3 only. In 
addition, we will take the dilaton $ to be a function only of these four- 
dimensional coordinates. The coordinates parameterize a. {d — 4)-torus 
with a reference radius tq. This radius will typically be chosen to be the 
string scale, or perhaps a few orders of magnitude larger. Similarly, 

the d-direction is compactified on a circle or an interval, either having a 
reference radius which, for simplicity, we also choose to be tq. For length 
scales much larger than this radius, all heavy modes will decouple and we 
will arrive at a four- dimensional effective theory describing the light modes. 
These light modes will necessarily include g'^^,, a'rp, o"^ and $. However, one 
should also make sure that there are no additional "stringy" light modes, 
such as winding modes around one of the internal circles. For the moment, 
let us ignore this important subtlety, returning to it at the end of this section. 
The associated action is most easily computed from (|Hn|) by Weyl rescaling 
the four- dimensional metric as 

V = ^^9iiu, (86) 

where 

= 6^*4 (87) 



and 



Then, making the field redefinitions 



= 2<I>1, dT = V2(c/ - 4)a^, <jR = V2a',,, (89) 

defining 

(90) 



1 VtVr 



where Vr and Vr are the volumes of the {d — 4)-torus and the (i-direction 
circle/interval respectively and dropping all higher derivative terms, we find 
that the low energy action is given by 

S=^J d'xV^C, (91) 



with 



C = R- ]^g^''d^<pd,<P - ]^g>"'d^aTd,aT - ^g'"' d^and^aR - 2e<^A. (92) 
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Note that this is exactly the four-dimensional theory introduced in ((T)) and 
if we set^ 

Q = /f • (93) 

Recall, however, that this result is predicated on the assumption that there 
are no additional stringy light modes. We must now check to see under what 
conditions this will be the case. 

To analyze this, we must be more careful in our discussion of decoupling. 
It is clear from (j85j) that, in the string frame (j83|) . the physical radii of the 
{d — 4)-torus and the d-direction circle are 

e"Tro, c^'rvq (94) 

respectively. From and (jH^ we see that, for the background geometries of 
interest in this paper, cr^ vanishes and, hence, the physical ((i— 4)-torus radius 
is time-independent and fixed at tq. On the other hand, it follows from (0]) 
and (j89|) that the physical radius of the (i-direction circle is time-dependent 
and given by 

rR(t) = |tanh(Qt)|ro. (95) 
For very early and late times, rfi{t) approaches tq. Note, however, that 

rnit) —.0, t-^ ±0. (96) 

Therefore, we must be concerned that stringy modes winding around the d- 
direction circle will become very light as t ^ ±0, thus substantially changing 
the low energy theory given in ()91|) and (j92|) . Can we estimate the time regime 
for which this effective theory is no longer valid? To do this, note that the 
winding modes around the ^-direction circle typically have a mass 

'^winding — ^ • 

(97) 

It then follows from ()95|) . and the fact that tq is of the order of the string 
scale or a few orders of magnitude larger, that the effective field theory (j^?T|l 

^In the simple supercritical string models we are considering, the cosmological constant 
A is of order the string scale, which is at odds with our assumption that Q is small 
compared to the momenta of interest. However, in the critical "cigar geometry" mentioned 
in the Introduction, it is possible to choose Q to be arbitrarily small. For this reason, we 
will simply ignore this issue in our analysis. 
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_____ — * 

and (|92p will be valid for momenta k satisfying 

r^tanh^(Qt) , , 

* « (ar 

Since we are interested in momenta satisfying^ 

Q2 < fc^ (99) 

we have to make sure that there is a window of momenta satisfying both (j^JHjl 
and (jnni)- Since we imagine Q to be much below the string scale, a'Q^ <^ 1, 
and for Tq = a', such a window exists as long as 

> a', (100) 

that is, as long as we stay more than a few string times away from the 
singularity. Note that this does not alter the discussion and conclusions 
associated with the theory given in and 0. To begin with, that theory 
is singular at the origin , so we could not really discuss the t —>■ ±0 regime, 
nor did we. Indeed, the stringy effects that descend from (j83p in the time 
region ()100p could conceivably "regulate" this singularity, making the theory 
well-defined at the origin. We will not explore this possibility in this paper. 
Henceforth, we will simply recall that the low energy theory (^, or, 
equivalently, (jHH), is not valid too close to the singularity. 

Subject to this caveat, we conclude that the four-dimensional cosmology 
described in this section descends from a class of string theories in + 1 
dimensions. We will now examine these string theories and their cosmology 
in detail, with the aim of eventually using them to compute the Bogolubov 
coefficients discussed above. 



3 String Theory Cosmology 

As we have just shown, the four-dimensional action (0), ((21) is the low en- 
ergy limit of the supercritical string theory (jH^]) dimensionally reduced to 

■^Massless modes with momenta below this scale correspond to growing modes and were 
discussed in section 5 of |15| . They lead to infrared divergences in string perturbation 
theory which will be ignored in this paper but would be interesting to understand better. 
We are not aware of any direct relation to the divergences in string perturbation theory 
mentioned at the beginning of section |3 
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four dimensions on a (c? — 4)-torus times a circle. The {d + l)-dimensional 
action (|83p is itself an effective theory, all higher a', string loop and non- 
perturbative effects being ignored. These effects are important, and we will 
comment on them in the appropriate places later in this paper. However, we 
will begin by studying the classical string action (jH^ and, specifically, the 
{d + l)-dimensional cosmological solutions of its equations of motion. One 
such solution, the "periodically identified generalized Milne solution" , or the 
"generalized Milne orbifold" for short, was first presented in jT3]. This turns 
out to be an exact solution of classical string theory, that is, it does not 
receive a' corrections We will see that, upon dimensional reduction to 
four dimensions, this solution includes the background spacetime studied in 
subsection 12 . II This makes the generalized Milne orbifold an excellent context 
to try and answer the questions raised in the previous section. 

3.1 The Background Spacetime 

Our starting point is the following solution to the equations of motion asso- 
ciated with the classical string action (j83|) . The metric and dilaton are found 
to be 

2 _ 1 dudv y^\^ 7^2 

$ = -^\og{l - uvf + %. (101) 

Here 

- oo < u, i; < oo (102) 
are global two-dimensional coordinates with the identification 

(m,i;) ~ (ue-'"«'-«,i;e2"^"°). (103) 

The parameter Q is related to the cosmological constant A by expression ()93|) . 
The first term in the metric describes the generalized Milne orbifold J3] . The 
remaining terms describe flat/torus-compactified {d — l)-dimensional space. 
The two-dimensional generalized Milne directions exhibit an intricate causal 
structure that is schematically represented in Fig. ^ 
In Regions I and II of Fig. ^ 

uv <0 (104) 
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VI 




II ' 



Figure 1: The generalized Milne orbifold [T3] . 

and the identification (jlU3|) is spacelike. These regions can also be described 
by the coordinates t, x defined by 



u = sinh(Qt)e-^^, v = - sinh(Qt) e^"" 

where 

a; ~ X + 27rro. 

In these coordinates, the metric and dilaton solutions become 

d-l 

1=1 

$ = - log cosh (Qt) + $0. 
For regions III and IV in Fig. ^ 

<uv < 1. 



(105) 
(106) 



(107) 



(108) 



The identification (jl03|) is now timelike. Thus, these regions contain closed 
timelike curves. In regions V and VI, one has 



1 < uv 



(109) 
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and the identification again becomes spacelike. 

Region I describes a circle that starts out at some fixed radius rg, collapses 
to zero size and then, in Region II, expands again to radius tq. The Big 
Crunch/Big Bang singularity occurs when 

uv = 0. (110) 

At this singularity, regions I and II touch the "whisker" [201 regions IV and 
III which, in turn, are connected to a second pair of asymptotic early and late 
time regions, VI and V respectively. In addition, the separation boundary 
between regions VI and IV and between regions V and III, defined by 

uv = 1, (111) 

is singular. Fig. [T] is intended to give a rough idea of the structure of the 
asymptotic regions, but is less precise about the structure near the singular- 
ities. For example, the metric in Region VI is actually such that the radius 
of the circle is increasing as one moves in towards the singularity. The string 
coupling Qs = exp($) grows as well. In drawing Fig. ^ we have implicitly 
performed a T-duality locally on that region, to bring it to a form more 
like that of Region I. Close to the singularity at uv = 0, where regions I, 
II, III and IV meet, the u, v piece of fllOip reduces to the Milne orbifold, a 
two-dimensional Minkowski space with the points related by fll03|) identified. 
See, for example, [211 EOl HH 121] for recent discussions. Deep into regions I 
and II, that is, as uv — *■ —oo, the dilaton becomes linear in time with the 
string coupling approaching zero. At uv = 1, the separation between regions 
VI and IV and between regions V and III, both the curvature and the string 
coupling diverge. 

A priori, one would expect stringy, higher derivative corrections to ()83p 
to become important near the singularities at -uf = and uv = 1. However, 
it was shown in ^H] that the background (jl(Jl|) defines an exact supercon- 
formal field theory with the correct central charge, c = 10 ^. That is, back- 
ground (jlOip is an exact solution of the string equations of motion to all 
orders in a'. 

The coordinates x^,x'^,x^ range from — oo to +oo. The coordinates 
x"',i = 4, . . . ,d — 1 are taken to be periodic, ~ x*^ + 27rro, so they de- 
scribe a (c? — 4)-torus with constant radii tq. These radii are chosen to be 

^The central charge deficit of tiie (u, v) directions is compensated by considering su- 
percritical string theory with a positive cosmological constant. 
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of the order of the string length or a few orders of magnitude larger. As 
mentioned previously, the coordinate x defined in (|1(J5|) is also compact al- 
though, in string frame, the corresponding circle has a time- dependent radius 
I tanh(Qt)|ro. The coordinate t will be the time coordinate of an observer 
living in regions I and II. 

We close this subsection by showing that in regions I and II of the {d + 1)- 
dimensional cosmological solution given in (jl(Jl|) . there exists an effective 
four-dimensional description of this background . This effective solution is 
valid everywhere except very close to the Big Crunch/Big Bang singularity, 
where it is expected to break down due to additional light stringy modes. 
It turns out that the dangerous modes are winding modes around the x 
direction, which, as we have previously discussed, indeed become light near 
t = 0. With this caveat in mind, we proceed to dimensionally reduce the 
metric and dilaton solutions in regions I and II from [d + l)-dimensions to 
four dimensions by compactifying them on the {d — 4)-torus times a circle. 
Recall that in regions I and II the metric and dilaton can be written as 
in (ll(J7|) . Using the definitions in subsection 12.41 we find that, at low energy, 
(ll(J7|) corresponds to 

ax = 0, 

an = v^log I tanh(Qt)|, 

= 2$o + log2-log|sinh(2gt)|, 

9f.u = e-V. (112) 

which is exactly the four- dimensional background given in Q, (HI), © and © 
of subsection 12.11 if we identify 

00 = $0 + ^ log 2. (113) 

It is crucial to note that ()l()7j) only describes regions I and II of the 
{d + l)-dimensional spacetime (jlOlj) . The full {d + l)-dimensional string 
background contains two additional asymptotic regions V and VI, as well as 
the intermediate regions III and IV. The four-dimensional spacetime ()112|) is 
a low energy description of regions I and II only. However, as we have pointed 
out before, this four-dimensional effective description breaks down near the 
Big Crunch/Big Bang singularity, which is exactly where regions I and II are 
connected to the other regions. To understand what happens near the Big 
Crunch/Big Bang singularity, as well as in regions III, IV, V and VI, it is 
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clearly necessary to use the full {d+ l)-dimensional string theory background 
(llUlj) . This is precisely what we will do throughout the remainder of this 
paper. 

3.2 Fluctuations 

In this subsection, we study fluctuations around the background (jlUl|) . The 
key ingredient in our discussion is that string theory allows one to determine 
globally defined wavefunctions, despite the singularities that prevent doing 
so in general relativity [T31 1201 • The underlying reason is that the spacetime 
(llUlj) corresponds to an orbifold of a coset conformal field theory.^ There 
is a well-defined procedure to determine at least the "untwisted" globally 
defined vertex operators in such theories. The zero mode parts of these 
vertex operators are the globally defined wavefunctions we are interested in. 

This procedure consists of two steps. In the first step, one determines 
the vertex operators of the coset conformal field theory, which in our case is 
PSL{2,JR)/U{1) at negative level.^ One describes the coset conformal field 
theory as a gauged Wess-Zumino-Witten (WZW) model. The vertex opera- 
tors of the ungauged PSL{2, JR) WZW model correspond to wavefunctions 
on the smooth PSL{2, IR) group manifold. These can be found in jSH]. The 
vertex operators of the SL{2,]R)/U{1) WZW model can then be viewed as 
those vertex operators of the ungauged model that are invariant under the 
[/(I) gauge group. See for example [1^1 1201 ESI- 

The second step is the standard string theory orbifold procedure jHTj . 
Roughly speaking, an orbifold is obtained from a covering space, in our case 
the PSL{2,'\R)/U{1) coset spacetime at negative level, by identifying points 
related by the action of a discrete group. Here, the group is Z with the 
action (jlOHj) . which turns a line into a circle (the circles visible in Fig. 
and causes the Big Crunch/Big Bang singularity at uv = 0. The untwisted 
vertex operators^ are those vertex operators on the covering space that are 

^For early work on applications of coset conformal field theories to cosmology, see for 
example [HSIISI- 

^PSL(2,]R)/U{1) at positive level corresponds to a two-dimensional black hole ge- 
ometry 35 . PSL{2,TR) /U{1) at negative level can be obtained from the black hole 
geometry by double Wick rotation 34 . It is described by the u, v piece of (|101|l , with- 
out the identification (|I03|I . In this spacetime, u = u = is a smooth point (the Big 
Crunch/Big Bang singularity only arises after the discrete identification (|1U3|) ). while 
there are singularities at uv — 1. 

^There are also twisted vertex operators, corresponding to strings winding around the 
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invariant under (|1(J3|) . From (jl(J6|) . we see that this amounts to momentum 
quantization in the x direction. In this paper, we will only be interested in 
zero momentum in the x direction since the x circle is part of the internal 
space. For the same reason, we will not be interested in winding modes 
except to note that they cause the effective four-dimensional description to 
break down near the Big Crunch/Big Bang. 

We would like to comment on the relationship of our theory to another 
model, which has the same Big Crunch/Big Bang singularity. This is the 
Milne orbifold (lR*-^'^^/Z) x IR^, where Z is generated by the boost trans- 
formation ()1U3|) on two-dimensional Minkowski space and denotes eight 
additional flat directions, some of which may be compactified. The associated 
metric is 



This spacetime is very similar to ()101|) near the Big Crunch/Big Bang sin- 
gularity u = V = 0, but differs significantly from it away from this point, 
in particular in the structure of the additional regions. The Milne orbifold 
consists of four cones touching at u = v = 0, that is, regions I, II, III and 
IV, but the radius of the circle of the cones grows to infinite size as one goes 
infinitely far away from the singularity. There are no regions V and VI. They 
can be thought of as having been pushed to infinity by focussing in on the 
manifold around u = v = 0. 

Wavefunctions on the Milne orbifold have been discussed from a string 
theory point of view in jSO] and have been used to compute string scattering 
amplitudes in [21]. There is no coset CFT involved in this spacetime and 
the untwisted wavefunctions can be easily obtained from the string orbifold 
procedure. That is, consider wavefunctions on Minkowski space and demand 
that they be invariant under the action (jl(J3|) of the orbifold group. Actually, 
such invariant Minkowski space wavefunctions either grow or decay in regions 
III and IV. Because these regions are non-compact in the Milne orbifold, 
one, often implicitly, further restricts to wavefunctions that decay in those 
regions. There will be no such additional restriction in the spacetime (I1U1|) . 
since regions III and IV do not extend to infinity there. As a consequence, 
more wavefunctions are necessary in our theory than one finds in the Milne 
orbifold. 

circles, see [TS] . 
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Global wavefunctions for the Milne orbifold with regions III and IV omit- 
ted were presented in [TT] based on a construction that does not refer to 
string theory, and applied in J2] to cosmology. The global wavefunctions 
agree with the restriction to regions I and II of the stringy wavefunctions 
discussed in the previous paragraph. Also, the vacuum state implicitly de- 
fined by the wavefunctions of [201 and more explicitly used in [21] corresponds 
to the vacuum state defined in ^Tj upon deleting regions III and IV of the 
string theory spacetime (or from the point of view of jTTj, upon adding those 
regions). It is of interest to see what happens to the wavefunctions of the 
generalized Milne orbifold upon taking the limit to the Milne orbifold. We 
analyze this in detail in subsection 2.4 and Appendix C. 

The fact that the wavefunctions that descend from a covering conformal 
field theory solution of string theory are globally defined on the associated 
orbifold is of fundamental importance for the results of this paper. For 
that reason, we outline in Appendix A, in more detail than discussed here, 
the procedure for constructing these global wavefuntions. In the remainder 
of this subsection, we will review and extend the results of on global 
wavefunctions and quantum vacuum states in the spacetime (jlUlj) . The 
global structure of the spacetime, in particular the presence of the additional 
asymptotic regions, will be shown to have important implications for the 
physics of the Big Crunch/Big Bang transition. 

3.2.1 Global Wavefunctions on the Generalized Milne Orbifold 

We restrict our discussion to those fluctuations that are relevant for the four- 
dimensional cosmology introduced in section |21 That is, the only non-zero 
momentum components are in the three non-compact space dimensions, and 
we ignore all winding and excited string modes. In particular, momentum 
and winding in the x direction are set to zero, 

p = w = (115) 

in the notation of 15 . The justification for this is that all the modes we 
ignore have masses on the order of the string scale, except very close to 
the Big Bang/Big Crunch singularity where winding modes around the x 
direction become light, as we have previously mentioned. 

In subsection \'2.'2[ we solved the wave equation (jHUj) for the fluctuations 
6T defined in (jSHJ. Recall that from the point of view of the four-dimensional 
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quantum field theory ((T]) and there were two classically independent re- 
gions, Region I and Region II. In each of these regions, we found the same 
two independent solutions of (jHUj) for a given momentum k. These solutions 
are STj:'^ in (jlH)) and 5T^~(= 5T^^*). Far from t = 0, one solution, 5T^^, 
reduces to a positive frequency wave, while its conjugate has negative fre- 
quency. However, as we have shown in subsection 13.11 this four-dimensional 
theory arises as the low energy limit of the generalized Milne orbifold solu- 
tion of the {d + l)-dimensional classical string action ()83p. The full string 
theory background (llOlj) has, in addition to regions I and II, the regions III, 
IV, V and VI discussed in detail above. Therefore, in string theory, one must 
solve the 6T fluctuation equation in each of these six regions. Recall that 
the generalized Milne orbifold arises from a covering PSL{2, IR) manifold by 
constructing the coset PSL{2,]R)/U{1) and then identifying points related 
by a Z group action. It follows from this structure that one can find solutions 
of the 6T fluctuation equation on the generalized Milne orbifold by first find- 
ing solutions of the fluctuation equation on PSL{2, IR) and then restricting 
to those solutions that are invariant under the action of both U{1) and Z. 
Such solutions, since they are globally defined on PSL{2,]R), remain glob- 
ally defined on the generalized Milne orbifold. That is, each such solution 
is defined in each of the six regions I, II, III, IV, V and VI. As explained 
in |15| 1^ . there are four independent wavefunctions of this type for each 
momentum k, which we denote by 

^++,fe(w^), IC_^j^{uv) and IC ^k{uv). (116) 

Note that, in addition to their k dependence, the argument of these functions 
is the coordinate uv. More precisely, they depend not only on uv but, also, 
on the region I, II, III, IV, V or VI. This latter dependence is suppressed in 
flllfij) . In regions I and II, it follows from p()5j) that 

uv = -smh^iQt). (117) 

Similarly, in regions V and VI we can write 

1 -Mt; = -sinh2(Qt). (118) 

Hence, in these regions the wavefunctions are dependent on t, as we would 
expect them to be. The four independent wavefunctions in ()116p can be 
written as 
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where is a fc-dependent normalization constant given in (j234j) and the 
expressions for each of the four functions Kj^j^ ^ in the six regions are given 
in Appendix B. That is, the independent wavefunctions in pi6|) are exphcitly 
known. The /C_|__|_ % are defined such that they are orthonormal, with norm 
squared ±1, with respect to the appropriate generahzation of the Klein- 
Gordon inner product (j43j) . The generalization is that S should be a global 
"Cauchy" surface, intersecting regions I and VI or 11 and V, which turn out 
to be equivalent. The Kj^^ g given in Appendix A have the same norm up to 
a sign, so, by dividing out a common factor, we can obtain the normalized 
/C^_i_ For example, /C^^ g coincides with 5T^ in Region 1 and vanishes in 
Region VI, so it indeed has unit Klein- Gordon norm. 

A key point of this paper is that the wavefunctions /C_,_^ ^ are globally 
defined on the generalized Milne orbifold. This follows from the fact that they 
descend from globally defined functions on the PSL{2, IR) covering space. 
The expressions for each of the four independent wavefunctions, in each of 
the six regions, are presented in Appendix B. It is, however, useful at this 
point to discuss several of these wavefunctions in more detail. First consider 
/C^^ In Region I, it is found to be 

jc^^rk = ^T^^^ (120) 

where 5^^^ is given in (PHj) . This function is purely positive frequency in the 
asymptotic region t <C —1/Q, while it diverges logarithmically near the Big 
Crunch singularity, t ^ —0. In Region II, one has 

k oc {-z)-^-'F{-j, -j; -2j- -), (121) 

where j and z are defined in pUj). In this region, is purely positive 

frequency asymptotically and diverges logarithmically near the Big Bang as 
t +0. In the intermediate Region III, 

g cx {uvyF{-j, -J- 1; 1 - — ), (122) 

which diverges logarithmically near the Big Crunch/Big Bang singularity at 
Mt; = 0, while it approaches a constant near the uv = 1 singularity that 
separates regions III and V. Similarly, in Region IV 

K^^ ,: oc {uv)-'-^F{j + 1, J + 1; 1; 1 - — ). (123) 
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This diverges logarithmically near uv — and approaches a constant near 
uv — 1. In Region V, /C^^ ^ is a mixture of positive and negative frequency 
components with equal amplitude in the asymptotic region uv ^ l/<5^, while 
it approaches a constant near uv = 1. The exact expression in Region V is 

}C++,k ^ + 1; 1; 1 - uv). (124) 

Finally, and importantly, in Region VI 

= 0. (125) 

The global behavior of /C^_ ^ can be obtained by replacing / VI, 
II V and uv ^ 1 — uv m. the above expressions for /C^^ Of particular 
importance for us is the form of /C^_ ^ in regions I and VI. We find that 

^+-,*? = (126) 

and 

in regions I and VI respectively. Similarly, the global structure of /C ^ is 

obtained from that of /C % by replacing / ^ J/, V ^ VI and "positive 
frequency" "negative frequency" . Specifically, we will use the fact that 



and 

^-^rt = (129) 

in regions II and V respectively. Finally, the behavior of /C_^ ^ is obtained 

from that of /C by replacing I <-> VI, II ^ V and uv -fr^ 1 — uv. In 

particular, we find that 

/C_^ g = (130) 

in Region II, whereas 

^-^,k - ST^^' (131) 

in Region V. 
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3.2.2 Quantization 



We would now like to quantize the scalar fluctuations bar on the generalized 
Milne orbifold. Generically, this can be done by expanding 



+a| ^K.\pv)e-^^-' + a\ ^¥il^^{uv)e-^^-') , (132) 



where {K.^ ^, K.^ is some appropriate pair of orthonormal wavef unctions. As 
compared to ()52p . we have twice as many functions in the expansion ()132|) . 
The reason is that there are twice as many independent global wavefunctions 
on the generalized Milne orbifold as there are independent wavefunctions in 
the effective four-dimensional theory. We now impose the canonical commu- 
tation relations 

K.A?' = - k'), a^,P^ = [a^ ay = 0. (133) 

and define the vacuum state |0) by 

ai|0) = aalO) = 0. (134) 

To proceed, we now must ask the question: what is a natural vacuum state 
to choose, or equivalently, what is a natural set of wavefunctions {/C^ IC^ g}? 
In ^3], two natural vacua were defined and shown to be inequivalent. The 
first vacuum state corresponds to the choice /C^ ^ = /C^_,_ K^^ = )C^_ ^. 
The associated fluctuation expansion is 

W^lCl^^{uv)e-^'^-^ + aY^lCl_^^^{uv)e~^^-^) . (135) 

Recall that at early times, /C^^ ^ is purely positive frequency in Region I 
and, from (|125p . vanishes in Region VI, and vice versa for the wavefunction 
/C , r. It follows that the vacuum state constructed from al and aY^ would 

H — ,k k k 

indeed be called empty by early time observers in regions I and VI. Therefore, 
we denote this state by |0)j„. The second vacuum state, specified by |0)o„i, 
corresponds to the choice /Cj^ ^ = ^*__ ^ and IC^j^ = ^ and the associated 
expansion 
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-ik-x 



Since, at late times, /C* 



j^[uv)e ^-^^ + a^'i^^_^^{uv)e '^'"j . (136) 
is purely positive frequency in Region II and, 



from p29|) . vanishes in Region V, and vice versa for /C_^ this state would 
similarly be called empty by late time observers in regions II and V. 
The relation between |0)j„ and \^)out can be determined from 
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\B* 
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V/c, 



— ,k 
-S 



(137) 



where the Klein-Gordon orthonormality of the wavefunctions implies that A, 
B and C satisfy the conditions 



\AY+\C\ 



\B\ 



1, 



AC* + A*C = 0. 



For the special case ()115|) we are considering, one finds 

A = -1, 
B = C = 



1 



^sinh(^) 



;i38) 



(139) 



Note that i? is a function of k and that B{k) = B{—k). Then, (I137|l simplifies 
to 
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which can be inverted to 
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(1C_ 
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It is straightforward to check these relations using the formulas in Appendix B. 
Inserting (|14H) into (jl35|) . we find the Bogolubov transformation 
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Note that the Bogolubov transformation (jl42j) mixes creation with anni- 
hilation operators. This imphes particle creation at late times ^H]- To see 
this, note that for each momentum k 

.„(0|af af|0).„ = \B\'= . ^ (143) 



siiih' (^) 



and 

,„(0|apa,^|0),„ =\B\'= ■ (144) 

smh -TTTT 



2Q 

The physical interpretation is that, at late times, the vacuum |0)j„ has 
1/ sinh^(7r£'g/2Q) particles with momentum k in Region I and the same num- 
ber of particles with momentum k in Region V for each value of k. Note that 
the number of particles created goes to zero^ in the limit that 

E-JQ oo. (147) 

We will use this result in subsection 13.41 

Let us now, within the context of the generalized Milne orbifold, compute 
the analogue of expression (|^. that is, the two-point correlation function 

in{0\SaT{uv,x)6crT{uv,x + r)\0)in. (148) 

Since this function is correlated with respect to the in-vacuum, it is simplest 
to expand the fluctuations as in (jl35p . Using this and (|133|) . we find that 



i{0\6aT{uv, x)5(Jt{uv, x + r)|0) 



in 



*ir(|K,,d^ + |K,_,^5iHWl). (149) 



2^2 j II V ++.'.1 ' ^-.f I 



^At least for the modes l)115|l we are considering. The generahzation of (|14H|I is [T^ 



2 cosh (7rm) 



(ff ) 

with 

n and w being integers labeUing momentum and winding in the x direction and E being 
the energy. 
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As in section 1, we will always set r = and consider 

.„(0|5aT(«t;,f)'|0).„ = J -0^ (|/C++,/ + 1^+-,/)- (150) 

Since our wavefunctions are globally defined, ()150|) is valid for {u, v) in any 
of the six regions of the orbifold. Of course, the explicit form of of /C+^ ^ and 

/C^ g change from region to region and, hence, so will the expression for the 

correlation function. Let us begin by calculating (jl5(J|) in Region I. It follows 
from (jT^ and (jT^ that 

,^{0\Sar{t,xy\OU = J (151) 

The agreement with ()5(i|l illustrates that the in-vacuum defined by pH5j) is 
indeed the correct vacuum in Region I. As previously, for momenta k"^ ^ 
in the distant past, the two-point function in Region I becomes 

where we have ignored a momentum-independent factor that can be removed 
by a field redefinition. The same expression for the correlation function and 
conclusions hold in Region VI. 

Having established this, we now calculate the object of real interest, 
namely, the two-point function ()150p in Region II. Using ()14H) and the ex- 
pressions for the wavefunctions given in the previous subsection, we find that 

in{0\5aT{uV,Xy\0)in 

= /(05((l + 2|5r)|5r^r + |5(fc)r5Tf + |5(fc)r5T^2).(153) 
As t ^ +00, this expression becomes 

,M6aT{uv,xr\0h^ = J 0^3j^{\ + ^)^ik,t), (154) 

where we have ignored a momentum-independent factor and 

A{k, t) = l + 2|5p + |5pe-2^^s* + |E| V'^s*. (155) 

This should be compared with the field theory result given in (f?^ . Clearly, 
it is impossible to find coefficients X and Y such that (f?^ reproduces ()155|) . 
This discrepancy has an interesting origin and interpretation, which will be 
the subject of subsection 13.31 
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3.2.3 A Family of Vacuum States 



So far, we have considered the two vacuum states of ^S]. One of them is 
empty in each of the early time regions I and VI, while the other is empty 
in the late time regions II and V. Using these vacua, we found that there is 
non-trivial particle creation in Region II, described by (jl43p . These particles 
influence the two-point function of the corresponding fleld, leading to the 
non-trivial factor of A{k,t) in ()155|) . Upon reflection, however, one is led to 
wonder whether the choice of vacua of [15j is really unique and, in particular, 
whether there exists alternative vacuum states for which the early time two- 
point function (jl5ip is less drastically altered upon going through the Big 
Crunch/Big Bang transition. 

Since we are interested in an observer who starts out in Region I and ends 
up in Region II after the Big Crunch/Big Bang transition, we will continue to 
impose the condition that any alternative in-vacuum state should be empty 
in Region I. That is, the two-point function in Region I should equal ()151|) . 
However, we will no longer impose a similar condition in the other in-region. 
Region VI, since our observer does not live there. We now construct a family 
of such generalized in- vacua and compute, in each case, the analogues of ()143p 
and ()154p . These correspond to quantities our observer could, in principle, 
measure in Region II. 

We continue to use the expansion (|132|) with 



which has positive frequency in Region I and vanishes in Region VI. However, 
we now allow a more general wavefunction for a-) which vanishes in Region I 
but does not necessarily have positive frequency in Region VI. The most 
general such wavefunction, for a given momentum, is an arbitrary normalized 
linear combination of /C+_ and }C*^__ specifled by 



l,k ~ 



(156) 



(157) 





For each 



l{k) ^ 0, 



(159) 
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there exists at least one 7 such that (jl58p is satisfied. Note that 7(A;) = 
1,7(A;) = corresponds to the natural in-vacuum |0)j„ defined using (jl35p . 
For simplicity, we will assume that 



(160) 



although more general choices would not change our conclusions. Henceforth, 
we will take 



/C 



2, A; 



77, fc 



(161) 



and use the expansion 
oai 



(27r)3/2 



k ++,k 
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ik-x ^ (jVJ^]Q* 
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ik-x 



■ (162) 



We can now repeat the computations of subsection 13.2.21 for the new 
in-vacuum states |0)^^ defined by ^H^ . (IT!?1) . (ITHH|l . (fTKTIl and (fTHKll . We 
continue to use the natural out- vacuum defined using p36|) . It follows from 
the expressions for the wavefunctions given in Appendix B and p57|l that 
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where i? is given in p39|) . Inserting this into ()162|) and comparing to ()136|) . 
we find the Bogolubov transformation 
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This can be inverted to give 
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Using (jl64j) . we can now compute the occupation numbers 

l7(fc)P 
(^) 



,(0|af af |0),, = |7i?P = . '7f/^JL (166) 



in Region II and 



,^(0|arta^|0),^ = |Sp + |7r (167) 



in Region V. For the reasons given at the beginning of this subsection, we are 
principally interested in the number of particles in Region II. We see from 
(jlfifjj) that this number can be made arbitrarily small by choosing 7 (A;) to be 
close to zero. Note, however, that we cannot set 7(/c) = 0, since this would 
be inconsistent with the constraint ()158p . 
Next we compute 

^;y{0\5aT{uV, x)^|0)^;y 

= / 1^ (|K,,,tP + . (168) 

In Region I, JC^ ^ vanishes. It follows that /C^- ^ also vanishes and, hence, 

(I168|) reduces to the familiar result (jl5ip . From (j225p in Appendix B, we see 
that /C_,__ j: is purely imaginary in Region II. In fact, this is also the case in 
the other regions that touch uv = 0. Therefore, /C^- g reduces to 7/C_,__ ^ in 
those regions. Using this, ()168p takes the following form in Region II, 

_fk_ 

(2^ 



(169) 



In the far future where t +00, and for k ^ Q , this expression becomes 

,,{0\5aAuv,xm,, = J (^l + ^^A{k,t), (170) 

where we have dropped a momentum-independent scale factor and 

A{k, t) = l + 2|75|2 + |75|2e-2^^£* + |75| V'^*^*. (171) 

In the limit '~f{k) —>■ 0, this approaches the early time result ()152|) . There- 
fore, for this limiting choice of in-vacuum, an observer in Region II would 
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say that the fluctuations in Region I went through the Big Crunch/Big Bang 
unchanged. For this special limiting case, the result (|171|) can be reformu- 
lated in the framework of subsection I2.31 It corresponds to the Bogolubov 
coefficients X and Y with the properties 

|X| = 1, Y = in the limit -f{k) 0. (172) 

However, this limiting case is the only one for which the framework of sub- 
section EISl can be made to reproduce ()171|) . This apparent discrepancy is 
the subject of the next subsection. 



3.3 Information Loss 

It is useful to use the Bogolubov transformation p65|) to explicitly write the 
relation between the in- vacua |0)^^ and the natural out- vacuum \0)out- We 
find that (remember p60p ) 

|0)^^ = iV^^exp jy" d^k {0afa^_!l: + Aa^a^t + fiai^^a^t)^ \0)out, (173) 

where 

e 

A 



2(7*(1 - 52) +7*)' 
7* - 7*52 

2(7*(1 - 52) +7*)' 



7*5 

To check this, use fll65|) to verify that |0)^;y deflned by ()173|) is indeed anni- 
hilated by and a^^ given that \0)out is annihilated by a^^ and a^. N^;y is 
a normalization constant, which we will determine for a special case. It is 
useful to keep in mind that B given in (jl39|) is purely imaginary. 

For the natural in- vacuum |0)j„ deflned in (jl35p . which corresponds to 
choosing 7(fc) = l,'j{k) = 0, the coefficients in p74|) become 





2(1 -fi2)' 
A = 

2(1-52)' 
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B 



(175) 



In this case, we find 
Nio = n (l - B'Ck))-'^" ^ exp [vj^^ log (l - B^k))-'^'^ , (176) 

with V the volume of space. For strictly infinity V, Niq vanishes, so we will 
assume V to be large but finite instead. Then ()176p can be derived as follows. 
Define the unitary matrix 

5 = exp {/ d'k (i?(af a!,^ + afa^_l - a^Ja^' - af^)) } . (177) 

Using the formula 

e-^5e^ = B + [B,A] + ^[[B, A], A] + ■ ■ ■ , (178) 

it is straightforward to verify that 

al = -Sal^S-\ ap = -5ap-^ (179) 

This implies that, up to a phase 

|0).„ = S\0)out. (180) 

The normalization factor p76|) can then be obtained by determining the 
coefficient of \0)out upon expanding the right hand side of (jlSOj) .^ 

As another special case, consider the limiting values 7(fc) —>■ 0, 7(A;) — *• 
oo. This was the limit that turned off particle creation in Region II, as we 
discussed at the end of subsection 13.2.31 This corresponds to 

^ ^ 0, 

/i ^ 0, (181) 



1 /2 

^In doing this, it is convenient to work with rescaled osciUators = (^{2tt)^/V) a^, 
which satisfy the canonical commutation relation [a^jjat] = 1. 
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and 

|0)ooo^iVo«,exp|y' ^^'^(^(aPa^t))} |0)out. (182) 

However, it turns out that in the hmit (jl81|) the normahzation factor iVg^oo 
vanishes, even for finite V. 

We have not computed the normahzation factor for generic values of 7 
and 7, although it could, in principle, be done in the same way. For instance, 
if 7 and 7 are both real and positive, the generalization of p77|l is 



S = expjy (f k p(a^i^^ a^j: + a- 

+/.(a,^ar,^-aPa^)}, (183) 



//f V] II yf II j 

k ~k k k k k ' 



with 



cosh(2/i) = 7 + 7, 
sinh(2/i) = —7, 

If 7 and 7 do not have the same phase, the exponent of 5* also involves cl^^O'^^^ 
and ar'^ a^l. 

k ~k 

The Hilbert space of states 7i we have been using can be viewed as the 
tensor product 

n = n^^0n^ (i85) 

of a Hilbert space H'' associated with a^^ and a^^^ and a Hilbert space 

associated with and a^. These two Hilbert spaces can be thought of as 
associated with the out-regions II and V respectively. Note that H could 
equally well be written as a tensor product of Hilbert spaces associated with 
in-regions I and VI. The out-vacuum can, accordingly, be expressed as 

\O)out = \O)ii0\O)v. (186) 



It is clear that each in- vacuum defined by ()173|) is a pure state in this tensor 
product Hilbert space. With a pure state, one can associate a trivial density 
matrix 

P77 = 10)7777(01- (187) 

When one is only interested in computing correlation functions in Region II, 
which would be the case for a physicist living in that region and trying to 
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predict the result of experiments done there, it is convenient (and equivalent) 
to use the density matrix obtained by tracing (|187|) over Ti.^ . This is given 
by 

plfi^ = Y.v{i\p^^\i)v, (188) 

i 

where {jz);/} is an orthonormal basis of TC^ . Note that |0)^;y G H = 7i^^ ^TiX , 
whereas |z)v G 7i^. Therefore, ()188p is indeed a density matrix in Region II, 
that is, an element of H^^* ® H^^ . Now, if p in ()173|) is nonzero (as it is, 
except in the limit 7 — » 0), |0)^;y is an entangled state in TC = H'' (8> . 
When one traces the density matrix of an entangled pure state over Ti^ , 
one obtains a non-trivial density matrix pSSj) in Region II. Therefore, an 
observer in Region II with no access to information about Region V finds 
himself in a mixed state, that is, a state with entropy. This entropy is 
called "entanglement entropy" and reflects the ignorance about the correlated 
Region V. In the limit 7 — 0, one again obtains a trivial density matrix 

PoL = |0)////(0|, (189) 

which is consistent with (jl72p . 

This explains the discrepancy noted at the end of subsection 13.2.31 In 
subsection 12.31 we assumed unitary evolution from Region I to Region II, 
that is, a pure state in Region I evolving into a pure (squeezed) state in 
Region II. This assumption generically turns out to be inconsistent with 
what we find in the full string theory background, which includes additional 
regions. We find that a pure in-vacuum state evolves to a pure state in 
the tensor product Hilbert space TC, but that this state generically contains 
correlations between regions II and V. As a consequence, in a description 
where Region V is ignored, such as the one appropriate for an observer in 
Region II, this state is a mixed state. Only in the limit 7 — > does the in- 
vacuum evolve to a pure state in Region II. As explained in subsection I3.2.3| 
the limit 7 — can, strictly speaking, not be reached. The more accurate 
statement is that the entropy of the state as described in Region II can be 
made arbitrarily small by taking 7 to be arbitrarily close to zero. 

The above discussion is a generalization of thermofield dynamics as re- 
cently applied to the eternal BTZ black hole, where the thermal state of a 
scalar field outside the black hole is obtained from an entangled state in a 
tensor product of two Hilbert spaces, the second Hilbert space being associ- 
ated with an additional asymptotic region behind the horizon of the black 
hole ^7j. Also see and references therein. 



47 



3.4 The Milne Orbifold Limit 



In the limit that 

uv -> 0, (190) 
the metric in (|1U1|) reduces to the Milne orbifold metric (|114|) 



1 ^ 

ds^ = —dudv + y{dx^f (191) 

Q 7=1 



subject to the identification ()103|) . 

(m,^;) ~ (Me-2-«^i;e2-«^). (192) 

In addition, it follows from (jlUlj) that, in this limit, 

$ — ^ $0- (193) 

Note that the theory we are explicitly discussing in this paper, (jlOll) . dif- 
fers in the number of compactified dimensions from the Milne orbifold ()114|) . 
Specifically, the generalized Milne orbifold has {d + l)-dimensions, whereas 
the Milne space has ten dimensions. How, then, can we claim that metric 
(jlOip reduces near the singularity to ()114j) ? The naive answer is that this 
limit does not involve the internal spatial dimensions which, therefore, are 
irrelevant. A more precise way of understanding the Milne limit is to rec- 
ognize, as we mention in the Introduction, that two fiat spatial directions 
of the generalized Milne orbifold can be replaced by a "cigar" geometry ^B] 
without changing any of the conclusion of this paper. In this case, we have 
a critical, ten-dimensional string theory where Q is a free parameter. This 
cigar geometry has a well-defined smooth limit to the Milne orbifold. Keep- 
ing this in mind, we will ignore the number of compact dimensions in this 
subsection. 

If in the limit (|19(jp we keep the ratio 

^ fixed, (194) 

then the fiuctuation equation (jHBj) on the generalized Milne orbifold reduces 
to 

5T" + -5T' + }?5T = 0. (195) 
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Were we to consider fields with non-vanishing momentum along the x circle 
as well, we would also be required to keep the product 

Qro fixed. (196) 

Expression (jl95|) is precisely the fiuctuation equation for a scalar mode in 
the Milne orbifold. The limit defined by (IMUl . (pi|l and takes the 

generalized Milne orbifold to the Milne orbifold. 

It is instructive to see what our global wavefunctions /C^_|_ g reduce to 
in this limit. In Appendix C, we show that /C_(_^ ^ reduces to the standard 
wavefunctions defining the adiabatic vacuum of the Milne orbifold, that is, 
superpositions of positive frequency plane waves in Minkowski space. They 
are Hankel functions in regions I and II and modified Bessel functions that 
decay asymptotically in regions III and IV, the expressions in the four re- 
gions being related by analytic continuation. Near the Big Crunch/Big Bang 
singularity, these wavefunctions diverge logarithmically. On the other hand, 
we find that 

(197) 

everywhere in the Milne limit. More precisely, it is always zero in Region I, 
is proportional to a finite Bessel function in Region II and corresponds to 
an asymptotically growing modified Bessel function in regions III and IV. 

Near the Big Crunch/Big Bang singularity, /C^ g approachs constants in all 

regions. However, its overall normalization factor vanishes, suppressing the 
wavefunction completely. It then follows from ()157|) that, in this limit, 

— 0. (198) 

Inserting these results into expression ()168|) for the scalar two-point function, 
we find that, in the Milne limit, for any choice of in- vacuum 

(^I^++,a;I'' (199) 

which is valid in all regions. This is the standard result for the two-point 
function in the adiabatic vacuum of the Milne orbifold. It is useful to evaluate 
this correlation function in regions I and II. It is straightforward to do this 
in Region I using (jl2(J|) . To evaluate (jl99j) in Region II, first note from (jl9(J|) 
and ()194j) that, in the Milne limit, 

E^/Q = ^JP-QVQ ^ oo (200) 
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and, hence, that the Bogolubov coefficient given in (|139|) satisfies 

B — ^ 0. (201) 

One can then compute the correlation function in Region II using (jl28|) and 
the Bogolubov transformation p63p . The result is that in both regions I and 
II of the Milne orbifold, 

,„(0|5aHt,f)^|0),. = / 0^,\STf\'. (202) 

It follows from our previous discussions that, in the distant past and dis- 
tant future, the spectrum is of the form of Minkowski fiuctuations plus a 
subdominant scale invariant contribution and, therefore, 

A(fc,t) = l. (203) 

That is, in the Milne orbifold, the spectrum is unchanged by passing through 
the Big Bang/Big Crunch singularity. Note that this quantum process can 
be described by the pure Region I and Region II four-dimensional effective 
theory with the Bogolubov coefficients 

|X| = 1, Y = 0. (204) 

Finally, we see that from ()143p and ()201|) that, in the Milne limit, 

,„(0|4'^^|0)„, — . 0. (205) 

That is, particle creation is turned off. We should mention that we have 
assumed that 7 is kept fixed during the limiting proceedure. It is clear that 
there would be particle creation if, instead, one kept \'yB\ fixed during this 
limit. 

3.5 A Note on Backreaction 

As we have emphasized at the end of the Introduction, the string theory 
background p01|) has been argued to be unstable to gravitational backre- 
action. Hence, reliable computations should take this backreaction into ac- 
count. This seems to be out of reach at present and we will not attempt to 
resolve the associated deep puzzles in string theory here. However, we would 
like to offer a few remarks from the point of view taken in this paper. 
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First, it is worth translating the issue into the four- dimensional language 
of section |21 Looking at the energy density (|12|) of the four- dimensional 
matter fields, one might be tempted to conclude that quantum mechanical 
particles should not lead to large backreaction. The reason is that the energy 
density of the classical matter fields, since it is dominated by scalar kinetic 
energy, scales like \t\~^ near the Big Bang/Big Crunch singularity at t = 0. 
The wavefunction of a particle diverges at most logarithmically near t = 0, 
so, according to (|^. its energy density cannot diverge more quickly than 
that of the background matter fields. It would seem, therefore, that a small 
fluctuation remains small compared to the background. So where is the large 
backreaction? 

The point is that the framework of section |21 is hard to use for com- 
putations near the bounce, since the geometry is singular there. The right 
framework to use is the one of sectional In these variables, there is not much 
energy in the classical matter fields near the bounce, where spacetime is lo- 
cally fiat with a slowly varying dilaton. The reflection of this in section |2] is 
that the sum of the energy densities in the matter and gravitational fields di- 
verges only like as one can see from (fT^ and (divided by nl). This 
sum is dominated by the potential energy of the four-dimensional dilaton 
which, in the string theory framework, arises from the cosmological constant. 
Therefore, the criterion for our computations to be unstable to backreaction 
is that the energy density of small fluctuations grows faster than l/\t\ near 
the singularity. Unfortunately, this is the case for generic fluctuations, whose 
wavefunctions diverge logarithmically at the bounce. Hence, we recover the 
familiar large backreaction problem. 

It is interesting to note that there is an important class of fluctuations 
that do not give rise to large backreaction |2I] . Of the modes considered in 
this paper, those described by /C^_ ^ are of this type, as can be seen using 
Appendix B. Their energy density actually vanishes at the Big Bang/Big 
Crunch singularity. 

The issue of backreaction is crucial from various points of view, string the- 
oretic as well as cosmological. A source of criticism of the simple orbifold-like 
cosmological singularities studied in the string theory literature is that they 
are unstable against even a single particle being added to the system before 

"'^'^This is not true for the K.^ ^ of the more generic string modes considered in |15j . 

For example, /C_| does not correspond to a chiral wavefunction near the singularity if 

it has non-zero momentum along the x circle. 
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the singularity [22] • Such a particle would change the conical singularity into 
a genuine curvature singularity. In our model, this is indeed the case for parti- 
cles with generic wavef unctions, as we have just discussed. Now suppose that, 
however remote the possibility, no such small fluctuation appeared and the 
Universe made it through the Crunch singularity. Then, we have shown that 
the Big Crunch/Big Bang creates a collection of particles with occupation 
numbers in Region II. One can wonder how these particles backreact 

on the geometry. First, will the energy density of these particles overwhelm 
the classical energy density very close to the Big Crunch/Big Bang singu- 
larity? The answer is no, as can be seen from ()168|) . The wavefunctions of 
the created particles are precisely proportional to /C^_ ^ near the singularity, 
so their energy density vanishes there. Second, it is tempting to speculate 
that, in some more realistic version of our model, this particle creation might 
play a role as a reheating mechanism. In fact, reheating by gravitational 
particle creation has been discussed before. See, for instance, jSH|- 
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Appendix A: Global Wavefunctions from String 
Theory 

In this Appendix, we will outline the procedure for constructing the globally 
defined wavefunctions on both the Milne orbifold and the generalized Milne 
orbifold discussed in this paper. 

^^Again, this would be different for the modes we ignored in this paper. In order to 
turn off their backreaction in Region II, one would have to work in a 7 ^ limit. 
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Milne Orbifold: We first review how globally defined wavefunctions are 
constructed on the Milne orbifold IR^'^/Z x IR^, an orbifold of Minkowski 
space. 

1. Solve the wave equation 

(V^ + m^)(p = (206) 

on the covering space IR^'^. A smooth basis of solutions is given by the 
plane waves 

|gip.Xg^(p+x-+p-x+)|2p+^- - / = m'}. (207) 

2. Choose an alternative basis of solutions to ()206p . consisting of continu- 
ous superpositions of the smooth solutions ()207|) . such that the orbifold 
generator i-^ e^'^^X'^ is diagonal. Such a basis is given by [3IT| (see 
also HU) 

0,^^. = e^P-^ r rf^e^(f+^-<=-'"+p-^+^™)e^'- (208) 

J — oo 

with I G IR. The orbifold generator acts by multiplication by e"^'^*'. 
The wavefunctions ()208|) are generically not smooth near the light-cone 
X~^X~ = 0, which divides the covering space in four regions. When 
studying string theory on Minkowski space, one could in principle use 
these singular wavefunctions but, since there exist smooth wavefunc- 
tions (j207p which can be written as continuous superpositions of (|208p . 
it is preferable to work with the latter. 

3. Restrict the wavefunctions ()208|) to the orbifold invariant ones, that is, 
to those with / G Z. Using coordinates t, x defined by = te^^/ y/2, 
in terms of which the orbifold generator acts a.s x ^ x + 27t, it is clear 
that / is the momentum in the x direction. Therefore, the condition 
Z G Z is the usual momentum quantization. Hence, we end up with 
wavefunctions ()208|) on the orbifold space. They are not smooth at 
the singularity X^X~ = 0, yet they are globally defined. Note that, 
because of the quantization of /, it is not possible to use a basis of 
smooth wavefunctions on the orbifold space. 

4. Introduce twisted sectors (winding modes). We will not discuss those 
in this paper; see |^ for a very recent discussion. 
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Generalized Milne Orbifold: We now review the prescription of [201 
IT^ for constructing global wavefunctions on the generalized Milne orbifold, 
which in string theory corresponds to a Z orbifold of a coset conformal field 
theory at negative level [T3] . 

PSL{2,]R)k<o . Tod-i 



U{1) 



/Z X IR'^"^ (209) 



1. Start with wavefunctions on AdS^, that is, the well-known eigenfunc- 
tions of the Laplacian on Ad S3 which coincides with the SL{2, JR) group 
manifold. SL{2, IR) is the group of 2 x 2 matrices g with unit determi- 
nant, 

g=(^ JjV ad~bc=l. (210) 



c d 

PSL(2, IR) is obtained from this by identifying g with —g. Wavefunc- 
tions on PSL{2, IR) should be consistent with this identification. The 
wavefunctions on a group manifold coincide with the matrix elements 
in the different representations in the group. Therefore, we can write 

^j,aA9) = {j,a\g\j,(3), (211) 

where j labels a PSL(2, IR) representation and a, (3 label states in the 
representation. In this paper, we will only consider the principal contin- 
uous representations, with j = — 1/2 + is and s real. They correspond 
to taking k"^ > Q"^. A choice of basis of wavefunctions now corre- 
sponds to a choice of states a,/3 in the representation j. By choosing 
an appropriate basis, one can choose a smooth basis of wavefunctions 
on PS'L(2,1R). Note that the group manifold PS'L(2,1R) admits a 
PSL{2,]R)l X PSL{2,]R)fi symmetry, the two factors acting on g by 
left and right multiplication respectively. It is clear from ()211|) that 
the action of PSL{2,'\R)i is determined by the state {j,a\, while the 
action of PS'L(2,1R)r is determined by The Laplacian on the 

PSL{2, IR) group manifold is the Casimir operator of PSL{2, IR) (act- 
ing either from the left or from the right). Its eigenvalue is j{j + 1). 

2. Choose a basis of wavefunctions on PS'L(2,1R) such that the U{1)l x 
U{l)n subgroup, with both factors generated by (T3, is diagonalized. 
This amounts to choosing the basis vectors of the representation j such 
that o"3 is diagonal. It turns out ^] that, for the representations we 
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are considering, there are two states for each eigenvalue m}'^ Thus, 
we find four independent wavefunctions for each j, m, rh 



One finds that the PSL{2, ]R) group manifold splits into six regions 
where these wavefunctions are smooth. Generically, the wavefunctions 
are not smooth at the boundaries separating these regions. However, 
they are still globally defined. When studying PSL{2, IR) by itself, one 
would usually prefer to work with a basis of smooth wavefunctions. 

3. The coset conformal field theory PSL{2,]R)k<:o/U{l), where k is the 
level of the PSL{2, IR) affine Lie algebra and U{1) acts a.s g ^ e°"' ge""^ , 
describes the generalized Milne spacetime ()101|) without the identifica- 
tion ()103p Its globally defined wavefunctions are obtained by 
restricting to those wavefunctions ()212|) that are invariant under 

that is, to those satisfying m = —rh [HHl ITH]. 

4. The Z orbifold group is another (discrete) subgroup of U{1)l x U{1)r. 
The globally defined orbifold wavefunctions are obtained by imposing 
a quantization condition on m (momentum quantization) and by intro- 
ducing twisted sectors (winding modes). We refer to [IS] for details. 
In the present paper, we are interested in modes with no momentum 
or winding along the Milne circle. 

Appendix B: More on Global Wavefunctions 

In this Appendix, we provide some technical details on the global wavefunc- 
tions g used in Sectional For simplicity of notation, we will suppress the 

subscript k. As compared to [131201, we restrict the discussion to the case 
of interest in this paper, as explained at the beginning of subsection 13.2.11 
In the notation of ^^1, this means that 



K±±{j, m, m; g) = {j, m, ±\g\j, m, ±). 



(212) 



A = /i = -J , 



(213) 



with j defined in and (UHl) as 




(214) 



12 



To be precise. 
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We start by recalling some properties of the hypergeometric function 
F{a, b; c; z) =2 Fi{a, b; c; z) Using the notation 

, , r(a + n) , , 

a „ = 215 
r(a) 

this function has the power series expansion 

F{a,b;c;z) = ±^^'-^, (216) 

n=0 ^■ 

which converges at least for |2;| < 1. It is useful to note that as z — >■ 0, 

F{a,b;c;z) — ^1 (217) 

for all values of the parameters a, b and c. The behavior for large \z\ can be 
obtained, for b — a not integer, using the transformation formula 

F{a,b-c-z) = pSp!^ ~ i i-zy^'Fia, 1 - c + a;l - b + a;-) 

r(o)r(c — a) z 

+fr^^r^(-^)"'^(^' ^-c + b;l-a + b;h (218) 
T{a)T{c — b) z 

which is valid for | arg(— z)| < tt. Equation ()218|) is singular for integer b — a. 
For a = 6, it is replaced by 

F(a a- c- z) = ^^'^ (-.)"» V Wn(l^c + a)^ 

na,a,c,z) r^^^p(^_^)l ^^ly 

X (log(-z) + 2^(n + 1) - ^(a + n) - ^(c - a - n)) ,(219) 

where 

^{z) = (220) 

Expression ()219p is valid for | arg(— 2;)| < vr, \z\ > 1 and c — a not integer. 
Other transformation formulas we will use are 

F(a, b; c; z) = (1 - z^F (^a, c-b;c; (221) 



and 



r(a)r(6) ^0 (n!)2 
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X (- log(l -z)+ 2^{n + 1) - *(a + n)- + n)) , (222) 

the latter of which is vahd for | arg(l — z)\ < n and |1 — -2| < 1- Also, recall 
that the Euler Beta function B{a, b) is defined by 

r>T6y ^^^^^ 

In terms of these functions, the global wavef unctions K±± are defined as 
follows. 



Region I: 

= ±-B{-j,-j){-uvyF l^-j,-j;-2j;-^^^ 

K._{uv) = ^B{j + 1,j + 1){-uv)-^-'f(j + 1,j + 1;2j + 2;—) 
2m \ uvy 

K+_{uv) = 

K^+{uv) = -B{-j,j + l)F{-j,j + l;l;uv) (224) 
m 

Region II: 

K++{uv) = ^B{j + l,j + l){-uv)-^-'F(j + hj + l;2j + 2;—) 
2m \ uvj 

K^4uv) = ^B{-j,-j){-uvyF (-j,-j;-2j;—) 
2m V uvJ 

K+_{uv) = -B{-j,j + l)F{-j,j + l-l;uv) 
m 

K^+{uv) = (225) 
Region III: 

K++{uv) = ^Bi-j,j + l){uvyF (-j,-r,l;l - —) 
2m \ uv / 



K^-{uv) = -^i?(-j,j + i)(H"'"'^(j + i,j + i;i;i- — 

1 / UV \ 

K+_{uv) = — s(-i,i + i)(i-H'^ 7 

ZTT? \ UV — IJ 

1 / UV \ 

K_+{uv) = B{-j,j + 1){1-uv)-^-'f(j + 1,j + 1;1; (226) 

2m \ UV — 1/ 

Region IV: 

K++{uv) = ^B{-j,j + l){uv)-^-'F (j + + —) 
2m \ UV/ 
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K^-{uv) = ^B{-j,j + l)iuvyF(-j,-r,l;l-—) 
2m \ uv J 

K+4uv) = —Bi-j,j + l){l-uv)-^-'F(j + l,j + l;l; 

Zm V uv — 1/ 

1 ■ / uv \ 

K.+ iuv) = B{-j,j + l){l-uvyF(-j,-r,l-l; (227) 

2m \ uv — 1/ 

Region V: 

m 

K^_{uv) = 

K+_{uv) = ^B{j + 1,j + 1){uv-1)-'-'f(j + 1,j + 1;2j + 2;-^ 
2m V 1 — uv 

K.+ {uv) = ^B{-j,-j){uv-iyF(-j,-j--2j--^) (228) 
2m \ 1 — uv/ 

Region VI: 

K++{uv) = 

K__iuv) = -5(-j,j + l)F(-j, j + 
m 

K+4uv) = ^Bi-j,-j)iuv-iyF(-j,-j;-2j; 



27ri V 1 — uv, 

K_+{uv) = ^B{j + 1,j + 1){uv-1)-^-^f(j + 1,j + 1-2j + 2--^) 
2m \ 1 — uv J 

(229) 

First, we use these explicit formulas to determine the behavior of some 
of these wavefunctions in the asymptotic early and late time regimes. In 
regions I and II, we can use (|1(J5|) to write 

uv = -sinh^iQt), (230) 

while in regions V and VI 

1 -Mf = -sinh^(gt). (231) 

Let us concentrate on the early time part of Region I (t ^ — 1/Q). It follows 
from (1221, and (1^ that as t -oo 

K++ - -j)2-2^e"2^«* = ^5(-J, -j)4-^e«*e-'^^^*. (232) 

58 



This shows that reduces to a positive frequency wave for early times 
in Region I. Moreover, expression (|232j) allows us to determine the factor M 
such that 

1C++ = (233) 

is normalized (see subsection 13.2. ip . The factor M is actually the same for 
all ir±±. Comparing (1^!?^ with (gTj) or (gHD, we find 

M = ^= . 234 

Using ()218|) . we can rewrite in Region I as 

K-+ = -B{-j,2j + l){-uvyF(-j,-j--2j-—) 
m \ uv J 

+ -B{j + 1, -2j - l){-uv)-^-^F (j + 1, J + 1; 2j + 2;—). (235) 

This shows that asymptotically in Region I, K |_ is a superposition of positive 

and negative frequency waves with equal amplitudes. 

Equivalently to ()235|) . we can rewrite Kj^_ in Region II as 



K+- = -B{-j,2j + l){-uvyF(-j,-j--2j- — 

HI \ uv / 

+ -B{j + 1, -2j - l){-uv)-^-'F (j + 1, J + 1; 2j + 2; — ) ,(236) 
m V uv/ 

which shows that asymptotically in Region II, is a superposition of 

positive and negative frequency waves with equal amplitudes. 

Next, we discuss the behavior near the Big Crunch/Big Bang singularity 
uv = 0. First consider In Region I, we can use (j219|) and (j224|) to write 

^++M = 7^f: ^~^¥n^^^" M"(-iog(-H+2^(^+i)-2^(-j+n)). 

n=0 v^v 

(237) 

We see that in Region I, diverges logarithmically near uv = 0. Similarly, 
we find in Region II 

K^^iuv) = — E r ^" M"(-log(-H+2^(^+l)-2^(j+l+^))- 

n=0 v""/ 

(238) 
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Applying (I^Hll and (1^^ to ((22^, we find that in Region III 



1 ^(-j)n(j + l) 



2m 



uv 



X (- log(Mt;) + 2^(n + 1) - ^(-j + n) - ^(j + 1 + n)). (239) 



The same equation, ()239p . holds in Region IV as well. Now consider 
near = 0. In Region I, it vanishes. The small uv behavior in Region II 
is manifest from (|216p and (j225p . In Region III, we can use (j221|) to rewrite 

(|22nD as 

K^^{uv) = ^^B{~j,j + J + 1; 1; uv). (240) 

The same formula, (j240j) . also holds in Region IV. Therefore, as uv ap- 
proaches zero, approaches a constant in regions III and IV, and twice 
that constant in Region II. 



Appendix C: Comparison with the Milne Orb- 
ifold 

In subsection 13.21 we gave a qualitative discussion of some differences and 
similarities between the pure and generalized Milne orbifolds. This compari- 
son can be made very concrete at the level of global wavef unctions, which we 
will do in this Appendix. Some of these observations were first made by the 
authors of As in the previous Appendix, we will suppress the subscript 
k for notational simplicity. 
In the limit that 

uv ^ with and Qvq fixed, (241) 

Q 

where vq is defined in p03p . the equation of motion (jHUj) reduces to that of 
a massless scalar in the Milne orbifold p95|) . We now take this limit of the 
solutions /C++ and /C_|__ and see what they correspond to in the pure Milne 
case. First, note that in the limit (j24H) . 

^ ^ +00. (242) 
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It then follows from (j^H) that 



1 



7 ^ h zoo. 



In this limit, the normalization factor 



becomes 



Using the fact that 

"^{z) ^ log(2;) for z — >■ oo with | arg(2;)| < vr, 
we find from ()237|) that in regions I and II 



(243) 



(244) 



(245) 



1 °° 1 



lot 



and 



1 °° 1 



log 



+ vri + 2^(n + 1) 



(246) 



■ni + 2^(n + 11 



respectively. Since in Region II uv = —Q t in the limit 



(247) 

the right hand 

r(2)/ 



side of fl247|) is precisely proportional to the Hankel function Hq (Et), that 
is, 

(^fc^), (248) 



while (I24(j|l is its analytic continuation to Region I with negative t via the 
lower half t plane JT]. Therefore, the limit ()24ip of /C++ in regions I and II 
is the wavefunction used in j^H HI] and jHO] to define the adiabatic vacuum 
inherited from Minkowski space. In Region III, using ()239p and ()243p . we 
find that 



K. 



1^1 ( Ejuv' 



++ 



E 



27ri — oin:: 



4g2 



log 



'Eluv' 



2^(n + 11 



(249) 



which is the expansion of the modified Bessel function Kq. 

1 



Til 





uvElX 


[i 





(250) 
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The function Kq decays asymptotically (that is, for large values of its argu- 
ment). This is exactly the analytic continuation of (|248|) to Region III and, 
additionally, to Region IV, via the lower half t plane. We conclude that in 
all four regions of pure Milne space I, II, III and IV, the limit ()24H) of /C++ 
exactly coincides with the wavefunctions used in |^ to define the vacuum 
inherited from Minkowski space. These wavefunctions coincide with those of 
[TT] when restricted to regions I and II. 

To determine the behavior of /Cj , note that in the limit 



f(-,., + l;l;»,)^i:|^(^] . (251) 

It follows that in Region II 

F{-j,j + l-l;uv)^Jo{E^t), (252) 
with Jo a Bessel function, whereas in regions III and IV 



F{-j,j + 1-1-UV) ^ lo 



uvEi 

k 



(253) 



with Iq a modified Bessel function which has the property that it grows 
asymptotically. For this reason, Iq it is usually not considered in discussions 
of the Milne orbifold, as we mentioned in subsection 13.21 In the generalized 
Milne orbifold, there is nothing wrong with growing behavior near uv = 
since regions III and IV do not extend to infinity. In the limit ()241|) . where 
regions III and IV are blown up to infinite size, the prefactor B{—j,j + 1) 
in (j225p and (j24(Jj) actually goes to zero. Therefore, although in regions III 
and IV is proportional to (j253p in the limit (j241|) . the proportionality 
factor is zero and K^^ vanishes, as it does in regions I and II. 
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